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I.  INTRODUCTION 

Over  the  last  few  decades  we  have  seen  an  explosion  in  the  advancement  of 
electronics  technology  which  has  resulted  in  smaller,  faster,  and  less  expensive  integrated 
circuits.  Integrated  circuits  are  the  building  blocks  of  today's  high-tech  systems  for 
commercial,  industrial,  and  defense  applications.  Our  increased  employment  of  space 
satellites  and  nuclear  power,  combined  with  shrinking  circuit  size,  dictates  the 
requirement  for  smaller,  more-sensitive  integrated  circuits  that  will  withstand  tougher 
radiation  environments.  Without  a  solid  understanding  of  the  radiation  environment, 
semiconductor  physics,  and  the  effect  of  radiation  on  semiconductor  devices,  our  ability 
to  operate  systems  which  affect  our  daily  lives,  such  as  communications  satellites,  is 
degraded. 

One  method  of  protecting  electronics  from  radiation  is  to  thicken  side  panels  or  to 
add  radiation  shielding  materials,  such  as  individual  shield  housings,  so  that  radiation 
penetration  is  decreased;  another  is  to  imbed  the  circuitry  within  other  components  to 
effect  radiation  shadowing.  These  options,  however,  may  not  always  be  possible  or 
practical.  For  most  systems,  such  as  spacecraft,  increasing  panel  thickness  and  adding 
individual  shield  housing  results  in  a  heavier,  more  costly  bus.  Additionally,  other 
factors  may  carry  more  weight  in  driving  a  system's  internal  configuration  than  radiation 
protection.  For  example,  in  a  geosynchronous  three-axis  stabilized  spacecraft,  the 
communications  equipment  is  mounted  on  the  north-  and  south-facing  bus  panels  for 
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thermal  isolation,  ruling  out  the  possibility  of  radiation  shadowing.   For  these  reasons,  it 
is  worthwhile  to  explore  a  third  option  of  utilizing  radiation  resistant  circuits. 

Electronic  circuits  may  be  built  from  components  which  are  radiation  hardened 
(rad  hard)  in  the  manufacturing  process.  Alternatively,  circuit  components  may  be 
arranged  into  novel  configurations,  or  have  inherent  characteristics,  which  increase  a 
circuit's  tolerance  (rad  tolerance)  to  radiation.  These  techniques  are  contrasted  in  Table  1 
[Ref  l:p.  2.4-1],  which  was  extracted  from  Harris  Semiconductor's  Space  Product  News. 
A  commercial  device  is,  as  might  be  surmised,  neither  rad  hard  nor  rad  tolerant. 

Rad  hard  components  are  guaranteed  by  the  manufacturer  to  survive,  i.e.,  remain 
within  data  sheet  limits,  when  exposed  to  a  given  radiation  level.  Dielectric  isolation, 
silicon  on  insulator  (SOI)  and  silicon  on  sapphire  (SOS)  technology,  and  higher  gate 
complexity  are  all  examples  of  radiation  hardness  design  techniques  used  by 
manufacturers  [Ref.  2:p  3.1.1-2].  While  rad  hard  techniques  are  the  most  reliable  of  the 
three,  the  complexity  of  the  manufacturing  process  may  prove  rad  hard  components  to  be 
unaffordable  for  some  applications,  particularly  in  an  era  of  budget  constraints. 

Rad  tolerant  devices,  on  the  other  hand,  provide  a  less  expensive  alternative  for 
protecting  system  circuits  in  a  radiation  environment.  Rad  tolerance  techniques  place  less 
emphasis  on  component  manufacturing  processes,  where  rad  hard  costs  are  incurred,  by 
employing  components  which  are  inherently  less  sensitive  to  radiation,  or  seeking 
innovative  ways  to  rearrange  the  commercial  components  within  existing  circuits  to 
minimize  radiation  sensitivity.  Naturally,  there  is  a  trade-off:  rad  tolerant  devices  are 


TABLE  1.1    GENERAL  CHARACTERISTICS  OF  RAD  HARD,  RAD 
TOLERANT,  AND  COMMERCIAL  INTEGRATED  CIRCUITS 


RAD  HARD 

RAD 
TOLERANT 

COMMERCIAL 

Designed  for  a  specific 

Hardness  offered 

Hardness  limited  by 

hardness  level 

as  a  by-product  of 

inherent  process  and 

the  design 

design,  customer  risk 

Wafer  lot  rad  test 

Sample  rad  test 

Customer  rad  test 

Guaranteed  to  remain  within 

Usually  tested  to 

Customer  rad  test  and 

data  sheet  limits 

functional  fail  only 

risk 

Total  Dose: 

20  krad  to  50  krad 

2  krad  to  10  krad 

>200kradto>lMrad 

(typical) 

(typical) 

Latchup 

Customer 

Customer  evaluation 

SOI  techniques  none 

evaluation  and  risk 

and  risk 

Bulk  rad  hard  techniques: 

guaranteed  extremely  low 

sensitivity 

less  reliable.  In  order  to  prove  radiation  tolerance  of  a  proposed  rad  tolerant  circuit, 
individual  lot-to-lot  testing  could  be  required  since  no  twin  components  are  precisely  the 
same.  Minute  manufacturing  differences  from  component  to  component,  which  might 
even  be  transparent  to  the  manufacturer,  may  lower  radiation  tolerance.  In  summary,  rad 
tolerant  devices  are  less  expensive,  albeit  less  reliable,  than  their  rad  hard  counterparts. 

The  operational  amplifier  (op  amp)  may  be  considered  to  be  the  basic  building 
block  of  a  wide  range  of  electronic  circuits.  An  op  amp  is  composed  of  numerous 
semiconductor  devices  which  may  be  radiation  hardened  to  collectively  protect  the  op 
amp  from  radiation  degradation.    Incidentally,  op  amps  composed  of  rad  hard  devices 
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inherently  have  large  parameter  fluctuations  as  a  result  of  the  rad  hard  manufacturing 
processes  [Ref.  3:p.  3].  The  response  of  op  amps  to  radiation  is  determined 
experimentally  for  a  single  application  when  exposed  to  a  particular  amount  of  radiation. 
If  the  circuit  is  later  modified,  further  testing  and  evaluation  are  required.  In  this  thesis 
single  op  amp  radiation  response  and  the  radiation  response  of  a  combination  of  op  amps 
which  replace  a  single  amplifier,  called  a  composite  operational  amplifier,  are 
investigated  and  compared. 

A  composite  op  amp  is  an  arrangement  of  two  or  more  op  amps  which  acts  as  a 
single  op  amp  but  with  improved  characteristics,  of  which  one  is  an  improvement  in 
radiation  sensitivity.  While  combining  single  op  amps  together  to  form  a  composite  op 
amp  will  not  increase  radiation  hardening  of  the  single  op  amp  circuit,  they  are  expected 
to  decrease  overall  active  circuit  sensitivity  to  radiation. 

Composite  op  amps  may  be  physically  realized  using  Bipolar  Junction  Transistors 
(BJTs),  Field-Effect  Transistors  (FETs),  or  Metal-Oxide  Semiconductor  Field-Effect 
Transistors  (MOSFETs  or  MOS).  Both  Bipolar  Junction  op  amps  and  Complimentary 
MOSFET  (CMOS)  op  amps  were  included  in  this  research,  which  simulates  and 
compares  the  radiation  response  of  single  and  of  composite  operational  amplifiers.  BJT 
components,  which  are  popular  due  to  their  fast  response  and  high  gains,  were  chosen  to 
allow  comparison  between  the  simulation  results  of  radiation  models  and  previous 
experimental  results  obtained  by  actually  irradiating  BJT  op  amps  using  a  linear 
accelerator  [Ref  11]. 


CMOS  op  amps,  on  the  other  hand,  are  popular  due  to  their  small  size,  high  input 
impedance,  easy  construction,  and  low  power  consumption.  The  drawback  to  CMOS 
devices  in  general  is  due  to  their  higher  sensitivity  to  small  amounts  of  ionizing  radiation. 
A  hybrid  circuit  utilizing  a  combination  of  BJT  and  CMOS  (BIMOS)  devices  which  is 
subject  to  ionizing  radiation  will  have  a  degradation  response  dominated  by  the  more 
sensitive  CMOS  devices.  For  this  reason,  it  was  important  to  include  CMOS  devices  in 
the  simulation. 

In  this  research,  single  and  composite  BJT  and  CMOS  op  amps  were  modeled  in  a 
finite  amplifier  configuration  using  the  PSPICE  circuit  simulation  program.  Total  dose 
radiation  effects  were  then  simulated  by  modifying  device  parameters  within  the  op  amp. 
The  outputs  of  both  composite  and  single  op  amps  were  then  compared  using  M ATLAB 
to  process  and  plot  the  extracted  data  from  the  PSPICE  output. 

Chapter  II  explores  natural  and  man-made  radiation  environments. 

Chapter  III  reviews  basic  semiconductor  physics. 

Chapter  IV  examines  radiation  effects  on  semiconductors. 

Chapter  V  reviews  op  amp  operation  and  composite  op  amp  theory. 

Chapter  VI  considers  radiation  effects  on  single  and  composite  op  amps. 

Chapter  VII  presents  the  models  used  and  steps  taken  in  the  simulation. 

Chapter  VIII  analyzes  and  displays  graphs  of  the  simulation  results. 

Chapter  IX  offers  recommendations  and  conclusions. 


II.  THE  RADIATION  ENVIRONMENT 

Radiation  may  be  traced  to  either  man-made  sources  or  natural  sources.  Natural 
sources  of  radiation  consist  of  the  sun  and  galactic  cosmic  rays,  while  man-made  sources 
of  radiation  of  interest  are  nuclear  reactors  and  nuclear  weapons,  all  of  which  are 
characterized  in  Table  2.1  [Ref  12:p.  4].  Although  they  are  quite  different,  it  is 
conceivable  to  have  both  sources  of  radiation  in  the  same  environment:     the  systems 


TABLE  2.1    SUMMARY  OF  MAJOR  RADIATION  DAMAGE  AREAS 


Application 


Damage  Effect 


Radiation" 


Device  Affected 


Displacement 


n  only 


Nuclear  weapons 


Ionization  • 


Transient         n  and  y 


Surface  n  and  y 


Bipolar  transistors 


All  semiconductor 
devices 


IGFETs 


Nuclear  reactor 


Displacement 


n  and  y 


Ionization  • 


Transient 


Bipolar  transistors 


No  problem 


Surface  y  only  Bipolar  transistors 

and  IGFETs 


Van  Allen  Belts 


Displacement 


eandp 


Ionization 


Transient 


Solar  cells  only 


No  problem 


Surface  ^andp  Bipolar  transistors 

and  IGFETs 


•n  —  neutrons 
e  —  electrons 


y  —  gamma  rays 
p  —  protons 


engineer  of  a  nuclear  powered  spacecraft  would  have  to  consider  both  sources  of 
radiation  when  configuring  equipment  panels  and  selecting  power  subsystem 
components.  For  this  reason,  it  is  necessary  to  understand  both  radiation  sources  and  the 
environments  they  create  in  order  to  predict,  and  ultimately  protect  against,  their  effects 
on  semiconductors.  In  this  chapter,  radiation-related  particles  and  processes  are  first 
catalogued  in  order  to  facilitate  the  ensuing  description  of  both  natural  and  man-made 
environments. 

A.    RADIATION-RELATED  PARTICLES  AND  PROCESSES 

1.    Charged  Particles 

Charged  particles  include  electrons,  protons,  beta  particles,  and  alpha  particles. 

Beta  particles  are  "fast  electrons"  emitted  from  decaying  nuclei  in  the  form 

n-+p++e~  and         p+^n  +  e+ 

and  so  they  may  be  further  defined  as  being  comprised  of  electrons  and  positrons,  also 

known  as  "antielectrons"  or  "positive  electrons."      Alpha  particles,  which  are  comprised 

of  helium  nuclei  (2  protons  and  2  neutrons),  are  much  larger  in  size  and  are  extremely 

heavy  compared  to  beta  particles.    They  are  the  result  of  nuclear  fusion  reactions.    Both 

alpha  particles  and  beta  particles  have  a  short  range  in  media,  such  as  air,  due  to  charge 

interactions.   This  is  particularly  true  for  the  relatively  sluggish  alpha  particles.    A  sheet 

of  aluminum,  such  as  a  satellite  panel,  or  even  the  next  collective  fifty  pages  of  this  thesis 

effectively  shield  most  alpha  particles.  [Ref  6:p.  143] 


2.    Photons 

Photons  are  small  packets  of  electromagnetic  energy  which  travel  at  the  speed  of 
light,  are  mass-less,  and  are  uncharged.  There  are  two  types  of  photons  of  primary 
interest:  X-rays  and  the  higher-energy  gamma  rays.  Both  are  products  of  nuclear  fission. 
X-rays  are  emitted  when  high-energy  charged  particles  incident  on  or  through  a  medium 
lose  kinetic  energy.  Gamma  rays,  which  are  highly  penetrating,  are  emitted  from  nuclei 
as  a  result  of  either  radioactive  decay  or  the  release  of  energy  after  nuclear  absorption  of 
a  neutron  (neutron  capture).  Photons  are  responsible  for  the  following  three  ionizing 
reactions: 

a.  Photoelectric  Effect 

When  a  relatively  low-energy  X-ray  passes  near  the  nucleus  of  an  atom,  it 
excites  the  atom  causing  it  to  shed  an  innermost  electron,  thereby  ionizing  the  atom  but 
losing  all  of  its  energy  in  the  process.  The  replacement  of  the  lost  electron  by  a  second 
electron  in  a  higher  shell  causes  the  emission  of  a  photon,  as  does  the  subsequent 
replacement  of  the  second  electron  by  a  third  in  an  even  higher  shell,  etc.,  known  as  the 
photoelectric  effect. 

b.  Compton  Effect 

The  Compton  effect  involves  the  collision  of  a  relatively  high-energy  X-ray 
with  an  electron,  resulting  in  ionization  and  only  a  loss  of  a  portion  of  the  photon's 
energy  (or  lowered  frequency)  instead  of  annihilation.  This  then  permits  the  X-ray,  and 
the  free  electron,  to  ionize  other  atoms  until  eventually  coming  to  rest. 


c    Pair  Production 

In  pair  production,  very  high-energy  photons,  or  gamma  rays,  are 
annihilated  when  passing  near  an  atomic  nucleus,  creating  an  electron  and  a  positron 
which  are  each  capable  of  ionizing  other  atoms. 

3.  Neutrons 

Neutrons  are  generated  by  both  nuclear  fission,  which  is  the  splitting  of  a  single 
element  into  multiple  elements,  and  nuclear  fusion,  which  is  the  reverse  of  fission,  and 
which  occurs  inside  the  Sun,  stars,  and  thermonuclear  reactions.  Both  reactions  release 
energy.  An  example  of  both  fusion  and  fission  is  the  fusion  of  a  neutron  and  a  235U 
nucleus  to  form  236U,  which  undergoes  fission  to  produce  fission  fragments  and  two 
prompt  neutrons:   [Ref.  18:p.  1 104] 

235^^236  U-*lAaXe+*Sr  +  2n 

4.  Ionization 

Ionization  damage  is  the  stripping  of  outermost  shell  electrons  of  an  atom  by  an 
incident  charged  particle.  In  semiconductors,  it  results  in  the  creation  of  free 
hole-electron  pairs,  such  as  in  the  three  photon  interactions  previously  discussed. 
Neutrons,  which  are  heavy  and  uncharged,  are  incapable  of  directly  causing  ionization; 
however,  they  may  induce  the  creation  of  other  ionizing  particles,  for  example,  neutron 
capture  by  a  nucleus  after  which  an  ionizing  gamma  ray  is  created.  Instead,  neutrons  are 
renown  for  displacing  entire  atoms  from  their  lattice  structures.  This  is  classified  as 
displacement  damage,  and  will  be  discussed  later. 


a.  Ionizing  Dose  Rate 

Ionizing  dose  rate  is  the  rate  at  which  ionization  is  accumulated  as  a  result 
of  the  deposition  of  charged  particles  and  photons  on  semiconductor  material. 

b.  Total  Dose 

Total  Dose  is  simply  the  integral  of  dose  rate:  it  is  the  accumulation  of 
ionization  a  circuit  receives,  measured  in  units  of  rads  of  silicon  (rads(Si)),  from  the 
absorption  of  charged  particles  and  photons  in  the  semiconductor  material.  One  rad  is 
one  hole-electron  pair  created  per  centimeter  squared  of  irradiated  material. 

5.    Displacement 

As  discussed  earlier,  displacement  damage  is  caused  by  the  displacement,  or 
dislodging,  of  an  atom  from  its  lattice  structure  by  relatively  heavy  particles.  Most  of  the 
damage  to  electronic  devices  in  this  environment  is  attributed  to  neutron  irradiation. 
Here,  the  relatively  heavy  neutrons  bombard  atoms  within  the  semiconductor  crystal, 
knocking  them  out  of  their  stable  positions.  These  displaced  atoms  may  either  take  up 
interstitial  positions  or  they  may  knock  other  atoms  from  their  positions  causing 
ionization  and,  ultimately,  permanent  defects  within  the  material  [Ref.  3:p.  161-2]. 
a.    Neutron  Fluence 

Neutron  fluence  is  the  rate  at  which  neutrons  are  incident  upon  a  material. 
Neutron  fluence  is  measured  in  neutrons  incident  per  centimeter  squared  of  irradiated 
material  (N/cm2). 
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6.  Single  Event  Upset  (SEU) 

Single  event  upset  is  the  upset  of  circuit  elements  due  to  a  single  particle  strike 
from  a  heavy  ion,  such  as  an  alpha  particle.  It  is  caused  by  the  deposition  of  charge 
within  a  critical  node  resulting  in  a  change-of-state,  such  as  changing  a  memory  cell  or  a 
logic  circuit  output  from  a  0  to  a  1  state,  or  vice  versa.  As  inferred,  it  is  a  serious 
problem  for  digital  circuits.  SEU  upset  is  not  included  in  this  thesis  for  two  reasons:  the 
first  is  that  SEU  is  caused  principally  by  heavy  ions  from  cosmic  rays,  and  is  not  a  major 
threat  in  orbits  lower  than  geosynchronous  since  the  earth's  magnetic  field  deflects  a 
majority  of  these  ions.  The  second  is  that  SEU  poses  a  threat  to  digital  circuits  whereas 
the  circuit  elements  and  the  proposed  technique  in  this  research  are  designed  for  analog 
circuits. 

7.  Annealing 

Annealing  is  the  ability  of  a  semiconductor  to  recrystallize  ("heal"  itself)  in  the 
presence  of  heat  or  current  after  having  experienced  radiation-induced  defects.  Because 
the  space  radiation  environment  is  a  continuous  process,  i.e.,  radiation  from  all  sources  is 
fairly  constant,  irradiated  devices  really  do  not  have  the  opportunity  to  anneal. 
Annealing  is  more  useful  in  the  context  of  a  pulse  of  radiation,  such  as  a  nuclear 
explosion,  after  which  an  irradiated  device  has  a  "rest"  period  during  which  it  may 
anneal. 
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B.    NATURAL  RADIATION 

Natural  radiation  exists  in  interplanetary  and  planetary  space.  Interplanetary 
radiation,  away  from  the  influence  of  the  earth's  planetary  magnetic  field,  is  caused  to  a 
small  degree  by  galactic  cosmic  rays  which  are  protons  and  electrons  from  stars,  and 
heavy  ions  from  sources  such  as  novas  and  supernovas.  Its  biggest  contributor  is  the 
solar  wind.  The  solar  wind  is  a  radial  stream  of  mainly  protons  and  electrons  emitted 
during  nuclear  fusion  within  the  Sun's  interior  whose  intensity  varies  with  solar  flare  and 
sun  spot  intensity  [Ref.  4:pp.  2.2.1-2]. 

Solar  flares  occur  randomly,  emitting  protons  and  a  few  alpha  particles  sparsely 
sprinkled  with  even  heavier  ions.  Radiation  from  a  solar  flare  reaches  the  earth  within  a 
few  hours,  peaks,  then  dissipates  within  a  week.  The  solar  wind,  however,  delivers 
residual  solar  flare  radiation  within  a  few  days  after  a  flare  through  striking  and 
compressing  the  Earth's  magnetosphere,  increasing  the  radiation  dose  at  lower  altitudes 
above  the  earth.  [Ref.  5:p.  II-8]. 

In  the  Earth's  planetary  space,  interplanetary  charged  particles  are  attracted  by  the 
Earth's  magnetic  field.  Those  particles  that  are  below  a  critical  energy  will  become 
trapped  and  gyrate  around  the  Earth's  magnetic  lines  of  force,  reflecting  between  mirror 
points  as  depicted  in  Figure  2.1  [Ref.  18:p.  695].  Collectively,  these  charged  particles 
form  bands  which  constitute  the  Earth's  Van  Allen  belts.  Figure  2.2  [Ref.  4:p.  2.2-4], 
which  is  an  extract  from  Harris  Semiconductor  manuals,  gives  a  rough  estimate  of  the 
radiation  exposure  rates  encountered  by  spacecraft  with  the  following  description: 
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Figure  2. 1   The  Earth's  Magnetic  Field 
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Figure  2.2    The  Spacecraft  Radiation  Environment 
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1 .  The  most  benign  environment  is  encountered  by  satellites  in  low  inclination,  low  earth  orbits 
(LEO)  shown  in  area  1 . 

2.  However,  low  earth  orbits  in  a  higher  polar  inclination  encounter  significantly  increased  dose 
rates  due  to  the  Van  Allen  belts  dipping  to  low  altitudes  and  a  much  higher  heavy-ion  flux  due  to 
the  lack  of  magnetic  deflection  of  the  ions  at  the  poles,  as  in  area  2. 

3.  In  area  3,  the  highest  dose  rates  of  up  to  1  Mega  rad  per  year  are  encountered.  This  area  is 
attributed  to  being  the  densest  region  of  the  Van  Allen  belts. 

4.  Spacecraft  in  a  geosynchronous  orbit  of  0°  inclination  and  35800  km  altitude  are  on  the  fringe  of 
the  Van  Allen  belts  and  experience  less  ionizing  radiation;  however,  at  this  altitude,  the  magnetic 
field  is  no  longer  strong  enough  to  act  as  a  shield  against  heavy  ions,  so  the  heavy  ion  flux 
increases.  [Ref.  4:p2.2-4] 

C.      MAN-MADE  RADIATION 

Unnatural  or  man-made  radiation  refers  to  nuclear  radiation  resulting  from  either 
nuclear  powered  systems  or  from  nuclear  weapons.  Radiation  from  nuclear  weapons  or 
nuclear  explosions  is  the  "worst-case"  in  that  it  is  characterized  by  infinite  dose  rates 
and/or  total  dose,  it  is  extremely  time-varying,  and  involves,  either  directly  or  indirectly, 
photons,  neutrons,  and  charged  particles.  For  illustrative  purposes,  the  following  section 
will  discuss  the  variety  of  particles  which  result  from  a  nuclear  explosion. 

Nuclear  radiation  is  released  both  at  the  instant  of  the  detonation  and  over  an 
extended  period  of  time.  Radiation  is  the  result  of  both  fission  and  fusion  reactions 
within  the  nuclear  device  and  is  released  both  at  the  instant  of  the  detonation  and  over  an 
extended  period  of  time.  Initial  radiation  consists  of  neutrons,  photons,  and  alpha  and 
beta  particles,  which  in  turn  cause  secondary  radiation  through  ionization  or  the  release 
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of  energy  after  capture.  Neutrons  are  generated  by  both  fission  and  fusion  detonation 
reactions  and  are  principally  responsible  for  producing  radioactive  materials  within  and 
without  the  explosion  confines.  Neutrons  travel  with  a  velocity  that  depends  on  their 
energy  and  interact  with  material  to  produce  additional  gamma  rays.  In  this  manner, 
neutrons  are  responsible  either  directly  or  indirectly  for  ionization  and  displacement 
damage.  [Ref3:pp.  142-3] 

Most  of  the  energy  of  an  explosion  is  emitted  as  X-rays,  which  are  the  principal 
agent  that  ionizes  the  air  to  create  what  is  known  as  a  fireball.  X-rays  are  readily 
absorbed  by  the  intervening  matter.  Gamma  rays,  on  the  other  hand,  are  very  fast  and 
highly  penetrating,  traveling  at  the  speed  of  light  regardless  of  their  energy.  Both  initial 
and  neutron-induced  gamma  rays  cause  ionization.  [Ref  3:pp.  142-3] 

Alpha  and  beta  particles,  as  discussed  earlier  in  this  chapter,  are  the  result  of  fusion 
reactions  and  decaying  nuclei,  respectively.  Due  to  their  short  range,  they  are  not 
considered  to  be  major  contributors  to  radiation  in  this  scenario.  [Ref.  3:pp.  142-143] 

D.    THESIS  RADIATION  ENVIRONMENT 

For  spacecraft  in  a  geosynchronous  orbit  and  below,  ionizing  radiation  trapped 
within  the  earth's  magnetic  field  poses  the  greatest  threat  to  electronics.  In  the  natural 
space  environment,  dose  rates  are  typically  low  (measured  in  rads(Si)  per  year)  and  are 
continuous  compared  to  other  radiation  sources,  such  as  nuclear  explosions,  where  a 
burst  or  pulse  of  radiation  is  experienced.  Radiation  incident  on  a  spacecraft  is  measured 
in  terms  of  accumulated  dose  at  the  end  of  the  spacecraft's  life  (EOL).    Spacecraft  are 
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designed  to  meet  mission  requirements  throughout  a  lifetime  total  dose  of  radiation.  For 
these  reasons,  total  dose  degradation  due  to  ionizing  radiation  (total  dose)  and  its  effects 
on  semiconductor  devices  are  simulated  and  are  the  main  subject  of  research  in  this 
thesis. 
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m.  SEMICONDUCTOR  DEVICES 

There  are  two  major  types  of  semiconductor  devices:  the  Bipolar  Junction  Transistor 
(BJT)  and  the  Field-Effect  Transistor  (FET).  The  majority  of  analog  circuits  are 
designed  using  BJTs  due  to  their  high  gain-bandwidth  product.  On  the  other  hand,  the 
use  of  FETs  is  minimal  in  analog  applications  because  of  its  relatively  small 
gain-bandwidth  product,  which  restricts  its  use  to  digital  circuits.  However,  the  FET 
offers  high  input  impedance,  easy  construction,  and  low  power  consumption.  Therefore, 
it  is  very  important  in  the  design  of  very  large  scale  integrated  circuits  (VLSI)  since  it 
may  be  made  quite  small  in  size,  thus  occupying  less  silicon  real  estate  on  an  Integrated 
Circuit  (IC)  chip.  This  chapter  first  discusses  the  BJT,  then  a  special  FET  circuit  called 
the  Complimentary  Metal  on  Oxide  Semiconductor  (CMOS)  FET. 

A.    BIPOLAR  JUNCTION  TRANSISTORS 

Within  the  BJT  there  are  three  semiconductor  regions:  the  emitter,  the  base,  and  the 
collector.  Each  region  consists  of  either  n-type  or  p-type  material.  N-type  material  is  a 
semiconductor  doped  with  atoms  that  can  donate  an  extra  electron,  thus  produces 
negative  carriers.  On  the  other  hand,  p-type  material  is  a  semiconductor  doped  with 
atoms  that  can  accept  an  extra  electron  which  creates  a  positively  charged  "hole,"  thus 
produces  positive  carriers.  For  an  npn  transistor,  both  the  emitter  and  the  collector  are 
made  of  n-type  material  and  the  base  is  made  of  p-type  material;  in  a  pnp  transistor,  the 
reverse  is  true.    Hence,  the  npn  and  pnp  BJT  may  be  respectively  modeled  as  either  two 
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back-to-back  or  two  front-to-front  diodes.  To  avoid  repetition  in  the  ensuing  discussion 
of  active  transistor  operation,  only  the  discussion  of  an  npn  transistor  will  be  considered, 
since  a  pnp  transistor  functions  similarly  except  the  current  is  reversed. 

A  transistor  is  biased  to  operate  in  the  active  linear  mode  if  it  is  to  operate  as  an 
amplifier.  An  example  of  an  npn  transistor  biased  as  such  is  shown  in  Figure  3.1  [Ref. 
5:p.  194],  where  current  flow  in  the  outside  circuit  represents  the  flow  of  positive  charge 
(holes).  Additionally,  Figure  3.2  [Ref.  5:p.  200]  illustrates  an  npn  BJT  cross  section. 


-      Vc»    + 


Figure  3.1  An  NPN  Transistor 


Figure  3.2  An  NPN  BJT  Cross-Section 
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1.  Emitter  Current 

Referring  to  Figure  3.1,  the  voltage  VBE  causes  the  p-type  base  to  be  higher  in 
potential  than  the  n-type  emitter,  which  forward  biases  the  base-emitter  junction.  This 
causes  a  current  to  flow  in  the  emitter  which  in  turn  is  the  sum  of  two  components:  holes 
transiting  from  the  base  to  the  emitter,  and  electrons  transiting  from  the  emitter  into  the 
base.  Due  to  a  higher  doping  of  n-type  material  in  the  emitter  than  p-type  material  in  the 
base,  the  hole  component  is  very  small  compared  to  the  emitter  component,  so  that  the 
emitter  current  refers  only  to  the  electrons  transiting  from  the  emitter  to  the  base.  These 
electrons  are  called  "minority  carriers,"  since  they  become  minorities  in  the  p-type  base 
region.  [Ref  5:pp.  194-195] 

2.  Base  Current 

Once  inside  the  base,  it  is  logical  to  assume  the  minority  carriers  recombine  with 
the  majority  carriers  within  the  base.  The  number  of  minority  carriers  that  recombine  is 
proportional  to  the  doping  level  of  the  p-type  base  material  and  inversely  proportional  to 
the  width  of  the  base.  In  actuality,  very  few  of  the  minority  carriers  recombine.  There 
are  two  reasons  for  this.  The  first,  which  was  alluded  to  earlier,  is  that  the  differences  in 
doping  between  the  emitter  and  the  base  lead  to  a  much  higher  population  of  electrons 
than  holes  in  the  vicinity.  Of  equal  importance,  the  width  of  the  base  is  purposefully 
very  thin.  So  the  cumulative  result  is  that  most  of  the  minority  carriers  that  transit  to  the 
base  are  able  to  continue  on  to  the  collector-base  depletion  region.  [Ref  5:pp.  194-195] 
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3.  Collector  Current 

Because  the  collector  is  biased  at  a  higher  potential  than  the  base,  i.e.,  V^  is 
positive  (reverse  biased),  those  minority  carriers  residing  in  the  base-collector  depletion 
region  are  swept  into  the  collector  as  collector  current. 

In  addition  to  the  above-mentioned  minority  carriers,  an  additional  parameter 
which  will  be  important  in  future  discussions  of  radiation  effects  is  reverse  currents. 
Thermally  generated  minority  carriers  produce  reverse  currents,  which  are  usually  small. 
However,  reverse  currents  in  the  collector-base  junction  are  actually  much  higher  than 
predicted,  and  are  dependent  on  both  V^  and  temperature.  [Ref  5:pp.  194-195] 

4.  Common  Emitter  Current  Gain 

Perhaps  the  most  important  transistor  parameter  is  the  common  emitter  current 
gain,  denoted  by  p.  While  (3  varies  from  transistor  to  transistor,  it  usually  ranges 
between  100-200  for  npn  transistors  and  a  little  less  for  pnp  type.  (3  is  the  ratio  of 
collector  current  to  base  current,  P=IC/IB  and,  as  such,  is  highly  dependent  on  the  relative 
doping  levels  of  the  collector  with  respect  to  the  base,  as  well  as  on  the  base  width. 
Since  J3  represents  current  gain,  it  is  desired  and  designed  to  be  of  high  value,  which  in 
retrospect  explains  why  the  base  is  thin  and  lightly  doped  in  comparison  to  both  the 
collector  and  emitter  [Ref  5:pp.  196-197]. 

B.    MOSFET  (MOS)  TRANSISTORS 

The  Field  Effect  Transistor  (FET)  is  also  referred  to  as  a  unipolar  transistor  since 
current  flow  is  due  only  to  the  motion  of  majority  carriers  in  contrast  to  BJTs,  whose 
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operation  depends  on  minority  carrier  flow  in  the  transistor  base  region.  FETs  also  differ 
from  BJTs  in  that  current  flow  is  parallel  to  the  surface,  whereas  in  a  BJT  it  is  mainly 
perpendicular.  [Ref  12:p.  141] 

As  the  name  implies,  a  Metal  Oxide  Semiconductor  (MOS)  transistor  has  an  oxide 
layer  inserted  to  insulate  the  metallic  base-like  terminal,  or  gate,  from  the 
semiconducting  material.  It  is  this  insulation  which  causes  the  current  in  the  gate 
terminal  to  be  extremely  small  and,  hence,  allows  the  overall  size  of  the  device  to  be 
small  [Ref  5.  p.  301].  Additionally,  MOS  devices  are  simple  to  construct,  so  that  large 
numbers  of  devices  may  be  produced  for  a  low  cost,  and  MOS  devices  are  also  known 
for  their  low  power  consumption.  A  CMOS  transistor  is  composed  of  two  MOS 
transistors:  one  which  consists  of  an  n-channel  imbedded  in  a  p-type  material,  called  an 
NMOS  transistor,  while  the  other  is  a  p-channel  imbedded  in  an  n-type  material,  called  a 
PMOS  transistor.  Below  is  a  description  of  NMOS  operation.  A  PMOS  transistor 
performs  similarly  to  an  NMOS,  only  the  current  flow  is  reversed. 

1.    NMOS  Transistor  Operation 

There  are  three  terminals  within  an  NMOS  transistor:  the  source,  the  gate  and  the 
drain.  Figure  3.3  [Ref.  5:p.  300]  provides  a  physical  schematic  of  an  NMOS  transistor. 
According  to  Sedra  and  Smith  [Ref.  5],  the  following  discussion  of  NMOS  operation  is 
divided  into  phases  which  are  a  function  of  the  magnitude  of  the  voltage  from  gate  to 
source,  vGS,  and  from  drain  to  source,  vDS. 
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Figure  3.3  An  NMOS  Transistor 

a.  Operation  Without  a  Gate  Voltage 

With  no  voltage  applied  to  the  gate,  vGS  two  back-to-back  diodes  exist 
between  the  drain  and  the  gate  terminals  which  will  not  conduct  current  in  the  presence 
of  a  potential  applied  between  drain  and  source,  vDS. 

b.  Operation  With  a  Small  Drain  Voltage 

With  the  source  grounded  and  a  small  positive  voltage  applied  to  the  gate, 
vGS,  holes  in  the  vicinity  of  the  gate  are  repelled  from  the  positive  gate  bias,  while 
electrons  are  attracted  and  fill  in  the  depletion  region  left  by  the  holes.  If  the  gate  voltage 
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vGS  is  of  a  sufficient  magnitude,  enough  electrons  accumulate  near  the  gate  to  form  a 
passage  between  drain  and  source,  thereby  creating  a  channel  for  current  to  flow.  The 
minimum  required  vGS  to  induce  current  flow  is  called  the  threshold  voltage,  denoted  by 
V, .  Vtis  typically  within  the  range  of  1  to  3  volts.  The  amount  of  current  that  flows  for 
a  given  vDS  is  controlled  by  the  magnitude  of  vGS  Hence,  the  greater  vGS  is,  the  more 
electrons  that  will  be  attracted  into  the  channel,  and  the  greater  the  channel  depth.  [Ref 
5:pp.  301-2] 

c.    Operation  With  a  Large  Drain  Voltage 

For  a  given  vGS  or  channel  depth,  an  increase  in  vDS  up  to  a  saturation 
voltage  of  vDS=vGS-Vt  results  in  a  linear  increase  in  the  current  flowing  through  the 
channel  as  depicted  in  Figure  3.4  [Ref.  5:p.  305];  however,  at  the  same  time,  the  channel 
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Figure  3.4  NMOS  Current 
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becomes  more  tapered  at  the  drain  end.  This  tapering  is  due  to  the  fact  that  vDS  appears 
along  the  length  of  the  channel,  varying  from  V=0  at  the  source  end  to  V=vDS  at  the  drain 
end.  Once  vDS  reaches  the  saturation  voltage,  the  drain  end  is  pinched  off  and  the 
transistor  is  said  to  have  entered  the  pinch-off  region,  where  vDS  no  longer  affects  current 
flow.  [Ref.  5:pp.  305-306] 

2.    CMOS  Operation 

CMOS  circuits  are  the  most  useful  of  all  the  analog  and  digital  integrated  circuit 
MOS  technologies.  A  CMOS  cross-section  is  illustrated  in  Figure  3.5  [Ref.  5:p.  307]. 
Note  that  while  the  PMOS  transistor  utilizes  the  n-type  body,  the  NMOS  draws  on  a  well 
of  p-type  material,  and  that  the  two  are  insulated  by  a  thick  region  of  oxide.     By 
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Figure  3.5  CMOS  Cross-Section 
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providing  both  transistor  types,  CMOS  transistors  allow  much  greater  flexibility  in 
circuit  design  compared  to  NMOS,  and  the  resulting  circuits  exhibit  improved 
performance  over  those  fabricated  with  NMOS  technology  [Ref  5:p.  A- 10]. 
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IV.  TOTAL  DOSE  EFFECTS  ON  SEMICONDUCTORS 

A.    DEGRADATION  IN  BIPOLAR  JUNCTION  TRANSISTORS 

While  new  BJT  devices  have  become  less  sensitive  to  displacement  damage  over  the 
past  ten  years,  total  dose-induced  gain  degradation  is  still  a  major  concern  for  bipolar 
transistors  [Ref  14:pp.  111-19-20].  This  is  particularly  true  if  BJT  devices  are  to  operate 
in  high-gain,  low-current  applications,  such  as  in  operational  amplifiers. 

Degradation  of  the  current  gain  parameter,  (3,  is  the  primary  effect  of  ionizing 
radiation  incident  on  bipolar  transistors.  Ionizing  radiation  penetrates  and  affects  both 
the  oxide  layer  and  the  substrate  beneath  it.  In  the  oxide  layer,  incident  radiation  creates 
inversion  layers  (positive  charge  on  an  n-substrate  surface  and  negative  charge  on  a 
p-substrate  surface)  along  the  exposed  transistor  surface  due  to  trapped 
ionization-induced  charge.  These  inversion  layers  increase  generation-recombination 
currents  and  create  space  charge  regions  along  the  oxide  layer.  This  phenomenon 
predominantly  affects  the  more-sensitive  base  region  due  to  its  relatively  thin  width  and 
low  doping  level.  The  overall  base  current  is  the  total  of  the  current  drawn  into  the 
base-collector  depletion  region  plus  the  current  drawn  into  the  inversion  layer  depletion 
region.  In  spite  of  increasing  base  current  and  shifts  in  the  collector-base  space  charge 
region,  the  much  larger  collector  current  (IC=(3IB)  is  impervious  to  these  relatively  small 
changes.  Increased  base  current  for  a  fixed  collector  current  reduces  the  bipolar 
transistor  current  gain  parameter  as  3=IC/IB,  and,  incidentally,  affects  the  transistor  most 
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when  the  transistor  has  a  high  gain  and  is  operated  at  low  current.  This  is  primarily  due 
to  the  fact  that  low  current  levels  are  surface  currents,  and  are  thus  affected  by  the 
inversion  layers  discussed  earlier.    Figure  4.1  [Ref.  3:p.  257]  depicts  the  relative  gain 

degradation  as  a  function  of  ionizing  dose  level  for  43  different  transistor  types. 
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Figure  4  1   Relative  Gain  Degradation  as  a  Function  of  Ionizing  Radiation  for  43 

Different  Transistor  Types 


While  these  effects  are  precipitating  from  radiation  incident  on  the  oxide  layer, 
ionizing  particles  which  weave  their  way  past  the  oxide  layer  are  able  to  shoot  into  the 
substrate,  creating  electron-hole  pairs  along  their  paths.  From  these  pairs,  the  highly 
mobile  electrons  contribute  to  the  minority  carrier  current  in  the  base.  Through 
recombination,  the  increased  base  current  shifts  the  base-collector  space  charge  region 
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further  into  the  collector,  increasing  the  base  width.    In  turn,  the  increased  base  width 
results  in  an  increase  in  base  transit  time  and  so  reduces  the  gain  bandwidth  product  ft. 
[Ref.  3:pp.  254-258] 

The  amount  of  base  current  generated  from  radiation-induced  inversion  layers  varies 
greatly  due  to  the  manufacturing  process  of  each  family  of  semiconductor  device,  which 
changes  the  device  structure  and  impurities  within  the  silicon  dioxide  layer.  A  particular 
semiconductor's  response  to  radiation  also  depends  on  the  bias  condition  under 
irradiation:  NPN  transistors  show  greater  gain  degradation  when  biased,  while  PNP  show 
greater  gain  degradation  when  they  are  unbiased.  Due  to  all  of  these  uncertainties,  a 
specific  correlation  between  (3  and  total  dose  is  unknown.    An  excerpt  from  Messenger 

[Ref.  3:pp.  255-256]  serves  to  legitimize  and  reiterate  the  above: 

Because  ionization  damage  ultimately  results  in  complex  surface  effects,  and  because  device 
structures  are  topologically  very  diverse,  there  is  at  present  no  theoretically  rooted  gam  degradation 
expression  as  a  function  of  ionizing  radiation  fluence  or  absorbed  dose,  as  there  is  for  bipolar  gain 
degradation  due  to  neutron  fluence. 

In  the  absence  of  equations,  there  are  some  general  guidelines  for  associating  (3  with 

total  dose  which  may  be  gleaned  from  a  few  references.    From  Messenger  and  Measel, 

[Ref  3  and  Ref  14],  a  reduction  in  (3  of  50%  results  from  an  irradiation  that  ranges 

between  .2  and  10M  rads(Si).  Figure  4.2  is  an  extrapolation  of  data  from  Figure  4.1,  and 

associates  a  reduction  in  (3  of  90%  to  a  total  dose  of  between  3M-50M  rads(Si),  and 

reductions  in  p  of  95%,  96%,  97%,  and  97.5%  to  total  doses  of  7M-70M,  7.5M-80M, 

8M-90M,  and  8.5M-100M  rads(Si),  respectively.    Referring  to  Figure  2.1,  10M  rads(Si) 
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Figure  4.1  Extension 
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Figure  4.2  Extension  of  Figure  4  1 


is  approximately  equal  to  the  amount  of  radiation  a  spacecraft  experiences  at  the  end  of  a 
ten-year  life,  while  orbiting  the  densest  region  of  the  Van  Allen  Belts,  which  is  at  an 
altitude  of  16,000  kms.  An  example  of  such  spacecraft  is  the  Global  Positioning  System 
(GPS)  satellites  which  orbit  the  Earth  at  approximately  this  altitude. 

With  all  of  this  in  mind,  it  is  important  to  emphasize  that  BJTs  are  much  more 
tolerant  to  ionizing  radiation  than  are  CMOS  devices.  Since  the  BJT  is  characterized  by 
small  base  widths  to  effect  high  efficiency  and  perpendicular  current  flow,  the 
above-mentioned  surface  effects  are  not  as  devastating  to  BJTs  as  they  are  to  CMOS 
devices,  which  depend  on  surface  interactions  for  their  operation. 

B.    DEGRADATION  IN  MOS  TRANSISTORS 

Shifts  in  MOS  transistor  threshold  voltages  due  to  silicon/silicon  dioxide  interface 
charge  trapping  is  the  dominant  mechanism  causing  device  degradation.   Figure  4.3  [Ref. 
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4:p.  2.3-2]  shows  a  schematic  representation  of  the  generation  and  trapping  of  positive 
charge  (holes)  within  the  oxide  layer  of  a  MOS  transistor.  Figure  4.3A  shows  the 
preirradiation  MOS  oxide  layer  cross  section,  while  Figure  4.3B  is  the  immediately 
post-irradiation  case.  Figure  4.3C  shows  a  reduction  in  electron-hole  pairs  due  to 
recombination,  where,  if  there  were  no  gate  bias,  all  of  the  pairs  would  recombine. 
Figure  4.3D  shows  the  mobile  electrons  migrating  to  the  positively  biased  gate  while  the 
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Figure  4.3  Generation  and  Trapping  of  Charge  in  an  NMOS  Transistor 
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much  slower  holes  migrate  away  from  the  gate  and  towards  the  channel  as  depicted  in 
Figure  4.3E.  Figure  4.3F  shows  hole  trapping  at  the  interface,  which  is  dependent  on 
interface  quality,  oxide  purity,  and  several  other  factors.  These  trapped  charges  cause 
transistor  thresholds  to  decrease,  such  that  N-channel  transistors  become  easier  to  turn 
on,  and  P-channel  transistors  become  harder  to  turn  on. 

Another  phenomenon  which  occurs  in  MOS  transistors  at  high  accumulated  doses  is 
interface  states.  Interface  states  are  generated  by  mechanisms  which  are  only  partly 
understood,  and  are  highly  dependent  on  fabrication  techniques.  Their  overall  effect,  as 
depicted  in  Figure  4.4  [Ref.  4:p.  2.3-3],  is  to  increase  N-channel  transistor  threshold 
voltages,  leaving  P-channel  thresholds  untouched.  Figure  4.4  shows  charge  buildup  and 
interface  state  effects  on  threshold  voltage  versus  total  dose  for  both  N-  and  P-channel 
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Figure  4.4  Charge  Trapping  Effects  on  MOS  Threshold  Voltages 
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MOS  transistors.  Due  to  the  varying  effect  between  P-  and  N-channel  transistors,  CMOS 
devices,  which  utilize  both  types,  have  complex  failure  modes  which  are  difficult  to 
predict.    For  these  reasons,  changes  to  CMOS  threshold  voltages  due  to  interface  traps 
would  be  very  difficult  to  analyze,  and  so  were  not  simulated  in  research  of  this  thesis. 
[Ref.  16:p.  1512] 
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V.  SINGLE  AND  COMPOSITE  OPERATIONAL  AMPLIFIERS 

An  operational  amplifier  (op  amp)  is  probably  the  most  widely  used  analog  circuit 
building  block  due  to  its  versatility  and  nearly  ideal  characteristics,  which  make  circuits 
utilizing  op-amps  easy  to  design.  There  are  generally  three  types  of  op  amps:  bipolar, 
CMOS,  and  a  combination  of  the  two  termed  BiCMOS.  Bipolar  op  amps  perform  better  in 
analog  circuits,  while  CMOS  op  amps  perform  better  in  digital  circuits. 

Made  up  of  a  large  number  of  transistors  and  resistors  and  a  single  stabilizing 
compensating  capacitor,  the  op  amp  is  a  complex  circuit  whose  performance  is  driven  by  the 
current  gain  within  each  transistor  ((3).  Figure  5.1  [Ref  5:p.700]  illustrates  a  741  internally 
compensated  op  amp  —  one  of  the  most  popular  op  amps  today,  which  is  currently  produced 
by  almost  every  manufacturer  of  analog  semiconductors.  As  such,  the  741  op  amp  is  the  type 
of  bipolar  op  amp  used  in  this  simulation. 

A.     OP  AMP  THEORY 

Driven  by  transistor  current  gain,  op  amp  performance  may  be  measured  in  terms  of 
three  parameters  which  will  be  discussed  in  turn:  3dB  frequency,  gain,  and  slew  rate.  The 
3 dB  frequency,  or  half-power  point,  is  important  in  that  it  is  a  measure  of  the  full  operational 
bandwidth  of  the  op  amp.  The  3dB  frequency  is,  in  turn,  determined  by  the  dominant  pole 
introduced  by  the  op  amp's  single  capacitor. 

Op  amp  gain  is  given  by  the  transfer  function  A=AocoL/(coL+s),  where  Ao  is  the  DC  gain, 
and  coL  is  the  3dB  frequency.   As  a  reminder,  coL  is  the  frequency  in  radians  per  second,  and  fL 
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Figure  5.1   741  Op  Amp 
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is  the  frequency  in  cycles  per  second  (Hz).  For  frequencies  »coL,  the  gain  reduces  to 
A=A0coL/s.  Another  parameter  is  the  gain  bandwidth  product  (ft).  The  f,  is  given  by  f=A0fL, 
which  is  a  constant  for  any  operational  amplifier  and  is  often  used  as  a  figure-of-merit. 
Figure  5.2  [Ref.  5:p.  724]  illustrates  A0,  3dB  frequency,  and  fr 
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Figure  5.2  Ideal  Op  Amp  Frequency  Response 

The  slew  rate  is  the  maximum  possible  rate-of-change  of  the  op  amp  output  voltage.  The 
slew  rate  depicts  the  ability  of  the  op  amp  output  to  follow  a  time-varying  input.  The  slew 
rate  of  a  compensated  op  amp  is  driven  by  the  time  needed  to  charge  or  discharge  the 
internal  compensating  capacitor.  Therefore,  the  slew  rate  is  relative  to  the  size  of  this 
capacitor,  which  is  needed  to  introduce  the  dominant  pole  in  the  op  amp  transfer  function  for 
stable  operation.  The  expected  time  response,  considering  only  the  dominant  pole,  is  a 
single-time  constant  reponse  given  by 


35 


Vout(t)=Vm(l-e,/tau) 

In  practice,  such  a  response  is  only  obtained  for  small  (less  than  lVpeak)  inputs.  For 
larger  inputs,  slew  rate  limiting  may  turn  a  sine  wave  into  a  triangular  wave  thus  introducing 
nonlinear  distortion.  The  maximum  slew  rate  of  an  op  amp  is  given  by  SR=Imax/C,  and, 
hence,  is  a  function  of  the  maximum  current  through  the  capacitor  in  the  second  stage  and  the 
magnitude  of  the  capacitance,  which  is  typically  small  to  effect  fast  slew  rates.  Op  amp  slew 
rate  is  determined  by  SR=fdVJdt\WiX,  and  is  usually  expressed  in  V/s. 

B.    COMPOSITE  OP  AMP  THEORY 

In  an  effort  to  extend  the  operational  bandwidth  and  improve  the  performance  of 
operational  amplifier  circuits,  Wasfy  B.  Mikhael  and  Sherif  Michael  developed  composite 
operational  amplifiers  (CNOAs,  where  the  N  denotes  the  number  of  operational  amplifiers)  in 
1981.  A  composite  op  amp  is  a  circuit  with  more  than  one  op  amp  that  permits  the  entire 
group  of  amplifiers  to  perform  as  a  single  amplifier  with  improved  characteristics.  [Ref 
ll:p.34] 

Initial  investigations  into  the  CNOA  behavior  have  been  discussed  in  the  literature  [Refs. 

6-11].  In  general,  the  procedure  used  for  developing  CNOAs  stemmed  from  creating  possible 

circuit  topologies  that  met  the  below-listed  performance  criteria: 
1. Of  the  136  possible  circuit  combinations  using  two  operational  amplifiers,  or  C20As,  the  denominator 
coefficients  should  satisfy  the  Routh-Horowitz  criterion  in  that  all  coefficients  be  of  the  same  sign- 
Additionally,  in  order  to  desensitize  the  C20A  with  respect  to  its  components,     no  numerator  or 
denominator  coefficient  should  be  realized  through  differences. 
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2.  The  terminals  of  the  C20A  should  resemble  as  closely  as  possible  those  of  the  single  op  amp. 

3 .  In  order  to  achieve  minimum  phase  shifts,  no  closed  loop  zeros  should  appear  in  the  right-half  of  the 
s-plane. 

4.  The  improved  frequency,  gam,  and  phase  performance  over  the  smgle  operational  amplifier  should  be 
enough  to  justify  the  mcreased  number  of  op  amps. 

From  the  criteria  listed  above,  the  C20A's  insensitivity  to  active  components,  passive 
components,  and  op  amps  makes  the  C20A  attractive  for  use  in  circuits  that  must  be  designed 
to  resist  radiation  degradation;  the  degradation  of  individual  circuit/op  amp  parameters  such 
as  gain,  slew  rate  and  3dB  frequency  will  have  less  of  an  effect  on  circuits  which  utilize 
C20As  than  on  those  which  use  single  op  amps.  Of  the  possible  136  C20A  circuit 
combinations,  the  four  listed  in  Figure  5.3  [Ref  6:p.  451]  proved  superior. 

Following  the  same  approach,  Composite  Multiple  Operational  Amplifiers  (CNOAs) 
were  developed  to  further  extend  the  operational  frequencies  at  the  expense  of  additional 
amplifiers.  C30As  were  obtained  by  starting  with  one  of  the  four  proposed  C20As  while 
replacing  one  of  its  single  op  amps  with  any  of  the  four  C20As.  C40As  were  developed 
similarly  with  each  of  the  C20A  single  op  amps  replaced  with  any  of  the  four  C20As,  or  by 
replacing  a  single  op  amp  of  the  C30As  by  a  C20A.  The  reader  interested  in  CNOAs  with 
N>2  are  referred  to  References  6  through  10,  as  only  C20As  were  analyzed  in  this  research. 
Figure  5.3  displays  the  circuits  of  C20A-1  through  C20A-4.  As  only  C20A-1  was  used  for 
comparison  in  this  thesis,  further  CNOA  discussion  will  be  focused  accordingly. 
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Figure  5.3  C20As 
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C.     COMPOSITE  VERSUS  SINGLE  OP  AMPS  IN  FINITE  GAIN  APPLICATIONS 

Finite  gain  amplifiers  are  used  in  instrumentation,  oscillators,  and  RC  active  filters. 
Ideally,  the  gain  of  these  amplifiers  should  be  constant  with  respect  to  frequency,  however, 
due  to  the  use  of  high  gain  op  amps  in  their  construction,  they  are  bandwidth  limited.  These 
op  amps  have  frequency  dependent  magnitude  and  phase  characteristics  which  can  cause  the 
circuit  response  to  change,  or  even  become  unstable.  One  way  to  overcome  this  problem  is  to 
improve  circuit  performance  by  the  use  of  actively  compensated  op  amps  that  have  superior 
characteristics,  such  as  CNOAs  [Ref  9:pp.  31-32]. 

Using  the  same  single-pole  representation  as  in  the  single  op  amp  gain  equation,  the  open 
loop  input-output  relationship  for  C20As  is 

VOI  =  VaAai(s)-VbAbi(s) 

where  i  denotes  the  type  of  C20A,  and  Aa  and  Ab  are  the  positive  and  negative  open 
loop  gains,  respectively.  It  may  be  shown  that,  for  the  C20A-1,  equation  1  corresponds  to 


Vol  =  VaA1{XA:A'){y  -  VbA-^± 


where  a  is  referred  to  as  the  compositing  resistor  ratio  of  C20A-1,  as  shown  above  in 
Figure  5.3.     If  both  single  op  amps  are  identical,  i.e.,  Al=A2=Ao5  and  the  C20A-1   is 

configured  for  a  single-ended  inverting  application,  then 
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Ao(C20A-i)=-2^  =  Ao(l+a)    for  (l+a)«Ac 


which  is  the  open-loop  dc  gain  for  a  C20A-1  inverting  amplifier.  The  C20A-1  will  have  a 
single  pole  roll-off  from  the  pole  at  wi/Ao  to  the  pole  at  wi/(l+ct).  Once  again,  it  can  be 
shown  that  the  transfer  function,  assuming  ideal  op  amps,  for  an  inverting  amplifier  is 

-k 


\+(s/wpOp)+(s2/wp2) 


where  »P=fig       and        Qp=  i/=gZ(l+<x) 

Low  sensitivity  to  circuit  parameters  is  guaranteed  since  none  of  the  numerator  or 
denominator  coefficients  is  realized  through  differences.  Additionally,  the  coefficients  in  the 
denominator  are  always  positive,  guaranteeing  stability.  Hence,  single  op  amps  with 
mismatched  gain  bandwidth  products  (fts)  within  practical  ranges  can  be  used  without 
appreciably  affecting  the  stability  or  sensitivity  of  the  finite  gain  realizations.  Simultaneously 
designing  for  a  maximally  flat  configuration  and  invoking  the  Routh-Horowitz  stability 
criterion  yields 


Op  =  -pr  and  a  + 1  > 


V  \+k 


where  a  =  J -^—  is  chosen,  and  k  is  the  overall  gain  of  the  circuit  utilizing  the  C20A-1 
[Ref6:pp.  449-455] 
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D.        BANDWIDTH     IMPROVEMENTS     USING     C20A-1     IN     FINITE     GAIN 
APPLICATIONS  AND  THE  IMPLICATIONS  FOR  RADIATION  TOLERANCE 

When  using  a  single  op  amp  for  a  finite  gain  amplifier,  the  bandwidth  shrinks  to 
approximately  (wi)/k.  Using  two  single  op  amps  cascaded  together  for  finite  gain 
amplification,  the  maximally  flat  3dB  bandwidth  is  obtained  when  each  amplifier  has  a  gain 
of 


A=fk 


in  order  to  realize  an  overall  gain  of  k.  The  resulting  bandwidth  shrinks  to 

A=  .66wi/  Jk 

The  C20A-ls  may  be  designed  to  shrink  to  only    w/Jk     for  the  maximally  flat  value  of 
Q=. 707.  Figure  5.4  [Ref  1  l:p.  48]  reveals  the  improved  bandwidth.  [Ref.  6:p.  455] 

The  large  increase  in  bandwidth  provided  by  composite  op  amps  could  be  important 

when  the  full  frequency  range  of  a  single  op  amp  is  required  during  its  operation  in  a 

damaging  radiation  environment.  The  composite  op  amp  would  have  to  undergo  severe 

damage  before  its  operating  frequency  range  is  reduced  to  that  of  a  single  op  amp. 

Combining  this  advantage  with  the  fact  mentioned  earlier  that  composite  op  amps  are  less 

sensitive  to  active  circuit  parameters,  it  is  clear  that  composite  op  amps  are  good  candidates 

for  use  in  radiation  environments. 
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Figure  5.4  Bandwidth  of  C20 As 
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VI.  RADIATION  EFFECTS  ON  OPERATIONAL  AMPLIFIERS 

Op  amps  are  relatively  sensitive  to  total  dose  effects.  Op  amp  gain,  3dB  frequency, 
and  slew  rate  are  some  of  the  parameters  affected.   [Ref.  19:  pp.  41-32-41-35] 

In  previous  research,  BJT  composite  and  single  operational  amplifiers  were 
irradiated  using  a  30  MeV  electron  beam  linear  accelerator  (LINAC),  where  in-situ  gain, 
slew  rate,  and  3dB  frequency  were  measured  to  ensure  irradiated  op  amps  were  not  given 
the  opportunity  to  anneal.  The  op  amps  tested  were  Harris  HS-5104RH  (rad  hard)  and 
HA-5104.  The  non  rad  hard  op  amps  are  a  general  purpose  op-amp  with  low  noise  and 
high  performance  (3  V/fis  slew  rate  and  8  MHz  gain  bandwidth  product),  and  they  were 
all  selected  from  the  same  lot.   [Ref.  1 1  :pp.  49-60] 

While  composite  op  amps  have  more  ideal  characteristics  than  single  op  amps  in 
general,  results  from  the  previous  LINAC  testing  (up  to  68M  rads(Si))  revealed  that 
composite  op  amps  also  degraded  at  a  slower  pace  (percentage-wise)  than  the  single  op 
amp  response  in  both  3dB  frequency  and  gain.  This  was  true  during  the  testing  of  both 
rad  hard  and  non  rad  hard  components.  The  difference  in  percentage  was  approximately 
5%  in  gain  and  15%  in  3dB  frequency  for  both  types  of  components.  Normalized  gain 
and  3dB  frequency  for  a  total  dose  of  6M  rads(Si)  for  non  rad  hard  op  amps,  and  18M 
and  68M  for  rad  hard  op  amps  are  displayed  in  Figures  6.1  through  6.3  [Ref.  ll:pp. 
72-73].  That  the  slew  rate  degradation  of  the  composite  op  amps  was  less  than  that  of 
the  single  op  amps  was  surprising,  as  it  was  predicted  that  the  degradation  would  be  the 
same.   [Ref.  1 1  :pp.  70-74] 
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Figure  6.1   Non  Rad  Hard  Op  Amp  Response  to  6M  rads(Si)  Using  a  LINAC 
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Figure  6.3   Rad  Hard  Op  Amp  Response  to  68M  rads(Si)  using  a  LINAC 
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During  the  LINAC  irradiation  test  of  single  and  composite  op  amps  in  the  previous 
research,  it  was  surprising  to  discover  that  the  slew  rate  of  the  composite  op  amp  was 
faster  than  of  the  single  amplifier  under  the  same  conditions.  It  was  also  surprising  to 
discover  that  the  rad  hard  devices,  composite  and  single,  survived  a  total  dose  of  over 
68M  rads(Si)  with  a  drop  of  only  5-10%  in  gain  and  a  50-60%  drop  in  3dB  frequency  (or 
bandwidth).  The  op  amp  performance  is  guaranteed  by  the  manufacturer  at  100k  rads. 

Incidentally,  Figure  6  1  is  a  prime  example  of  how  unique  device  response  to  radiation 
can  be:  even  the  LINAC-irradiated  twin  op  amps  on  the  same  chip  did  not  degrade 
similarly 
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Vn.  SIMULATION  PROCEDURE 

To  simulate  the  effect  of  radiation  degradation  on  op  amp  parameters,  the  following 
technique  is  introduced.  The  main  concept  in  this  technique  is  to  generate 
correspondences  between  1 )  transistor  parameter  variations  and  total  dose  radiation,  and 
2)  transistor  parameter  variations  and  circuit  response,  so  that  a  combination  of  the  two, 
in  turn,  would  correlate  the  relationship  between  circuit  response  and  total  dose  radiation. 
The  correspondence  between  transistor  parameter  variations  and  total  dose  radiation  was 
already  described  in  Chapter  IV,  specifically  in  Figures  4.2  and  4.4.  A  correspondence 
between  transistor  parameter  variations  and  circuit  response  was  obtained  through 
iterating  the  same  inverting  amplifier  configuration  response  several  times  using  different 
transistor  parameter  values,  for  which  the  procedure  is  discussed  in  this  chapter. 

A.    OP  AMP  MODELS 

The  first  step  was  to  choose  both  BJT  and  CMOS  op  amp  models  for  use  in  both  the 
single  and  composite  op  amp  configurations.  Two  models  were  selected  from  Sedra  and 
Smith  [Ref.  5]  due  to  their  simplicity,  popularity,  and  availability  of  their  devices' 
PSPICE  parameters.  The  BJT  and  CMOS  op  amp  circuit  layouts  used  in  this  simulation 
are  displayed  respectively  in  Figure  7.1,  which  was  adapted  from  Figure  5.1,  and  Figure 
7.2  [Ref.  5:p.  C-6],  followed  by  their  associated  PSPICE  programs. 
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Figure  7. 1   BJT  Op  Amp  Circuit 
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♦741    OP  AMP    PSPICE    CIRCUIT 


.SUBCKT  OPAMP  5,14,19 

VCC 

1  0  15 

VEE 

26 

0  - 

-15 

Ql 

12 

5 

8 

MOD2 

Q2 

12 

14 

7 

MOD2 

Q3 

10 

6 

8 

MODI 

Q4 

15 

6 

7 

MODI 

Q5 

10 

13 

11 

MOD2 

Q6 

15 

13 

16 

MOD2 

Q7 

1 

10 

13 

MOD2 

Q8 

12 

12 

1 

MODI 

Q9 

6 

12 

1 

MODI 

Q10 

6 

3 

4 

MOD2 

Qll 

3 

3 

26 

MOD2 

Q12 

2 

2 

1 

MODI 

Q13 

17 

2 

1 

MOD  3 

Q14 

1 

17 

22 

MOD5  3 

Q15    17    22    19    MOD2 

Q16    1       15    9      MOD2 

Q17    18    9      25    MOD2 

Q18    17    20    21    MOD2 

Q19    17    17    20    MOD2 

Q2  0    26    21    23    MOD6    3 

Q21    24    23    19    MODI 

Q22    15    24    26   MOD2 

Q23    26    18    21    MODI 

Q24    24    24    26   MOD2 

Q25    18    2       1      MOD4 

Rl       11    26    IK 

R2       16    26    IK 

R3       13    26    50K 

R4       4      26    5K 

R5      2       3       39K 

R6      22    19    27 

R7       19    23    27 

R8       25    26    100 

R9       9      26    50K 

R10    20    21    40K 

Rll    24    26    50K 

CI       15    18    30P 

.MODEL   MODI    PNP(BF=50    IS=10F    CJE=.1F    CJC=1.05P    CJS=5 

.MODEL   MOD2    NPN(BF=200    IS=10F   CJE=.65F   CJC=.36P    CJS= 

.MODEL   MOD3    PNP(BF=50    IS=2.5F    CJE=.1F    CJC=.3P    C JS=4 . 

.MODEL   MOD4    PNP(BF=50    IS=7.5F   CJE=.1F   CJC=.9P    C JS=4 . 

.MODEL   MOD5   NPN(BF=200    IS=10F    CJE=2.8F   CJC=1.55P    CJS 

.MODEL   MOD6    PNP(BF=50    IS=10F    CJE=4 . 05F    CJC=2.8P) 

.ENDS    OPAMP 
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Figure  7.2  CMOS  Op  Amp 


*CMOS    OP   AMPPSPICE    CIRCUIT 

.SUBCKT    CMOSOPAMP    2,1,7 
VDD      10    0    DC      5 
0    DC    -5 
9    DC    25U 


VSS 

IREF 

Ml 

M2 

M3 

M4 

M5 

M6 

M7 

M8 

CC 

RC 

CL         7 
.MODEL    1 
.MODEL    2 


13  10  2  L=8U  W=120U 
2  3  10  2  L=8U  W=120U 
4  9   9   1  L=10U  W=50U 

4  9   9   1  L=10U  W=50U 
6  10  10  2  L=10U  W=150U 

5  9   9   1  L=10U  W=100U 

6  10  10  2  L=10U  W=150U 

6  10  10  2  L=10U  W=150U 
8  10PF 

7  10K 

0  10PF 

NMOS  (LEVEL=2  VTO=l  KP=20U  LAMBDA=0.04) 

PMOS  (LEVEL=2  VTO=-l  KP=10U  LAMBDA=0.04 


.ENDS  CMOSOPAMP 


51 


B.  OP  AMPS  IN  THE  INVERTING  AMPLIFIER  CONFIGURATION 

The  inverting  amplifier  configuration  was  chosen  because  it  is  simple  and  provided  a 
comparison  with  previous  research  (Sage).  Since  radiation  degradation  in  transistors  is 
dependent  on  bias,  it  was  important  to  have  roughly  the  same  input  voltage  as  Sage's 
circuit  to  produce  a  useful  comparison.  Sage's  input  voltage  was  4mV  Vp  with  an 
overall  gain  of  200.  A  circuit  with  an  input  voltage  of  lOmV  Vp  AC  and  inverting 
amplifier  gain  of  100  for  an  output  gain  of  IV  was  chosen  in  this  thesis. 

C.  RADIATION  SIMULATION 

1.    Gain  and  3dB  Degradation 

In  simulating  ionizing  radiation  in  both  BJT  and  CMOS  transistors,  the  current 
gam  (PSPICE  code  BF)  and  threshold  voltages  (PSPICE  code  VTO)  were  changed, 
respectively,  within  each  transistor  of  each  circuit.  (3  was  varied  from  0%  to  97.5% 
degradation  in  both  NPN  and  PNP  transistors,  corresponding  to  respective  drops  from 
200  to  6,  and  100  to  3.  Each  BJT  op  amp  circuit  simulation  was  executed  seven  times 
with  decreases  in  P  of  0%  50%,  90%,  95%,  96%,  97%,  and  97.5%.  In  order  to  preserve 
data  points  for  the  nonlinear  portion  of  the  degradation  curves,  a  large  discontinuity 
between  P  degradation  of  0%  and  90%  was  allowed  since  little  change  in  op  amp 
response  was  noted  between  these  data  points.  The  value  of  97.5%  was  selected  as  the 
end  point  since  P  degradation  beyond  this  point  destabilized  the  circuit,  for  both  single 
and  composite  op  amps. 
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In  the  CMOS  op  amps,  Vt  was  altered  from  IV  to  -IV,  and  from  -IV  to  -3V  for 
NMOS  and  PMOS  transistors,  respectively.  Since  the  CMOS  components  are  more 
sensitive  to  ionizing  radiation,  with  smaller  changes  in  parameters,  more  regularly  spaced 
Vt  degradation  data  points  were  selected  than  for  the  BIT  transistors.  Degradation  points 
of  Vt,  Vt-.25V,Vt-.5V,  Vt-lV,  Vt- 1.25  V,  Vt-1.5V,  and  Vt-2V  were  selected,  after  which 
the  circuit  for  both  single  and  composite  op  amps  became  unstable. 

An  AC  analysis  was  conducted  in  PSPICE  for  each  inverting  amplifier  circuit  which 
corresponded  to  a  particular  degradation  value  of  (3  and  of  Vt.  The  programs  for  single 
and  composite  BJT  and  CMOS  op  amps  are  given  below.  After  running  each  circuit  for 
each  data  point,  the  voltage  output  was  extracted  into  two  sets  of  data,  with  frequency  as 
the  first  column,  and  output  voltage  (in  dB)  at  node  3  in  the  second  column.  Keeping  in 
mind  that  CMOS  data  processing  paralleled  that  of  the  BJT  data,  programs  for  loading 
and  processing  only  BJT  data  are  provided  below.  Data  for  all  circuits  was  gathered  in  a 
MATLAB  file  similar  to  the  program  entitled  MATLAB  LOADING  PROGRAM  and 
arranged  in  a  matrix  where  rows  one  through  fourteen  contained  single  op  amp  circuit 
data  (frequency  and  voltage)  for  each  of  seven  3  degradation  points,  and  rows  fifteen 
through  twenty-eight  contained  associated  composite  op  amp  data. 

This  large  data  matrix  was  transferred  to  a  second  MATLAB  file,  such  as  the 
program  entitled  MATLAB  GAIN  AND  3dB  PROCESSING  PROGRAM  below,  which 
processed  the  matrix  to  obtain  gain  and  3dB  frequencies,  and  then  saved  the  result.  To 
obtain  the  inverting  amplifier  finite  gain  for  each  frequency  interval  in  dB,  the  input 
voltage  of  lOmV  was  converted  to  dB,  then  subtracted  from  the  output  voltage.  The  3dB 
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*PSPICE    BJT    SINGLE    OP   AMP     (741)     INVERTING 
* CONFIGURATION    SHOWING    BETA   CHANGES 

.SDBCKT   OPAMP 


*0%    DEGRADATION 

.MODEL   MODI    PNP (BF=50    IS=10F    CJE=.1F   CJC=1.05P    CJS=5 . IP) 
.MODEL  MOD2    NPN(BF=200    IS=10F   CJE=.65F   CJC=.36P    CJS=3 . 2P) 
.MODEL  MOD3    PNP (BF=50    IS=2 . 5F   CJE=.1F   CJC=.3P    CJS=4 . 8P) 
.MODEL   MOD4    PNP (BF=50    IS=7 . 5F   CJE=.1F    CJC=.9P    CJS=4 . 8P) 
.MODEL   MOD5    NPN(BF=200    IS=10F   CJE=2 . 8F   CJC=1 . 55P    CJS=7 . 8P) 
.MODEL   MOD6    PNP (BF=50    IS=10F    CJE=4.05F   CJC=2 . 8P) 
%50%    DEGRADATION 

.MODEL  MODI    PNP(BF=25    IS=10F   CJE=.1F    CJC=1 . 05P    CJS=5 . IP) 
.MODEL  MOD2    NPN (BF=100    IS=10F   CJE=.65F    CJC=.36P    CJS=3 . 2P) 
.MODEL  MOD3    PNP(BF=25    IS=2.5F   CJE=.1F   CJC=.3P    CJS=4 . 8P) 
.MODEL  MOD4    PNP(BF=25    IS=7.5F   CJE=.1F    CJC=.9P    CJS=4 . 8P) 
.MODEL  MOD5   NPN(BF=100    IS=10F   CJE=2 . 8F   CJC=1.55P   CJS=7.8P) 
.MODEL  MOD6    PNP (BF=25    IS=10F   CJE=4.05F   CJC=2 . 8P) 
*90%   DEGRADATION 

.MODEL  MODI    PNP(BF=5    IS=10F    C JE= . IF   CJC=1 . 05P    CJS=5.1P) 
.MODEL  MOD2    NPN(BF=20    IS=10F    CJE=.65F   CJC=.36P    CJS=3 . 2P) 
.MODEL  MOD3    PNP (BF=5    IS=2 . 5F    CJE=.1F   C JC= . 3P    CJS=4.8P) 
.MODEL  MOD4    PNP (BF=5    IS=7 . 5F    CJE=.1F   CJC= . 9P    CJS=4.8P) 
.MODEL  MOD5    NPN(BF=20    IS=10F   CJE=2 . 8F   CJC=1 . 55P   CJS=7 . 8P) 
.MODEL  MOD6    PNP (BF=5    IS=10F   CJE=4 . 05F   CJC=2 . 8P) 
*95%    DEGRADATION 

.MODEL  MODI    PNP (BF=2 . 5    IS=10F    C JE= . IF   CJC=1.05P    CJS=5 . IP) 
.MODEL  MOD2    NPN(BF=10    IS=10F    CJE=.65F   CJC=.36P    CJS=3 . 2P) 
.MODEL  MOD3    PNP (BF=2 . 5    IS=2 . 5F   CJE=.1F    CJC=.3P    CJS=4 . 8P) 
.MODEL  MOD4    PNP (BF=2 . 5    IS=7.5F    CJE=.1F   CJC=.9P    CJS=4.8P) 
.MODEL  MOD5    NPN(BF=10    IS=10F    CJE=2 . 8F   CJC=1 . 55P    CJS=7 . 8P) 
.MODEL  MOD6    PNP (BF=2 . 5    IS=10F    CJE=4 . 05F   CJC=2 . 8P) 
*96%    DEGRADATION 

.MODEL  MODI    PNP (BF=2    IS=10F    CJE=.1F    CJC=1 . 05P    CJS=5 . IP) 
.MODEL  MOD2    NPN (BF=8    IS=10F   CJE=.65F   CJC=.36P    CJS=3 . 2P) 
.MODEL  MOD3    PNP (BF=2    IS=2.5F    CJE=.1F    CJC= . 3P    CJS=4.8P) 
.MODEL  MOD4    PNP (BF=2    IS=7.5F    C JE= . IF   CJC= . 9P    CJS=4 . 8P) 
.MODEL  MOD5    NPN (BF=8    IS=10F    CJE=2 . 8F   CJC=1 . 55P    CJS=7 . 8P) 
.MODEL  MOD6    PNP (BF=2    IS=10F   CJE=4.05F    CJC=2 . 8P) 
*97*    DEGRADATION 

.MODEL  MODI    PNP (BF=1 . 5    IS=10F    CJE= . IF   CJC=1 . 05P    CJS=5 . IP) 
-MODEL  MOD2    NPN(BF=6    IS=10F    CJE=.65F   CJC=.36P    CJS=3.2P) 
.MODEL  MOD3    PNP (BF=1 . 5    IS=2 . 5F    CJE=.1F   CJC=.3P    CJS=4 . 8P) 
.MODEL  MOD4    PNP (BF=1 , 5    IS=7.5F   CJE= . IF   C JC= . 9P    CJS=4 . 8P) 
.MODEL  MOD5    NPN(BF=6    IS=10F   CJE=2 . 8F   CJC=1 . 55P    CJS=7 . 8P) 
.MODEL   MOD6    PNP (BF=1 . 5    IS=10F    CJE=4.05F   CJC=2 . 8P) 
ENDS    OPAMP 
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*PSPICE    BJT    SINGLE    OP  AMP    (741)     INVERTING 
♦CONFIGURATION   SHOWING   BETA   CHANGES  CONTD. 


♦INVERTING      AMPLIFIER    CIRCUIT 
XI         0    2    3         OPAMP 

VIN      4    0    AC    .01  *gain   analysis   AC    input 

VIN    4    0    PULSE (0     .01    10U    IN    IN    . 13M    .26M)*slew   rate    input 
Rin      2    4    le3 
Rout   2    3    le5 

*PLOT   ANALYSIS 
.AC   DEC    40    le3HZ    le6HZ  *gain  AC   analysis 

.PRINT   AC    VDB(3,0) 
.PLOT      AC    VDB(3,0) 

. TRAN    5U    . 3M  *slew   rate   analysis 

.PRINT    TRAN   V(3)    V(4) 
.PLOT      TRAN   V (3)      (0,-1)    V(4)      (0,.01) 

.END 


♦COMPOSITE  (C20A1)  BIPOLAR  OP  AMP  INVERTING  CONFIGURATION 

. SUBCKT  OPAMP  *refer   to   single  hj t   op   amp  program 


.ENDS  OPAMP 

*C20A1  CIRCUIT 
. SUBCKT  C20A1  1,2,3 

XI  2  4  5  OPAMP 

X2  1  5  3  OPAMP 

Rl  4  1  IK 

R2  5  4  6. IK 
.ENDS  C20A1 

♦INVERTING   AMPLIFIER  CIRCUIT 
XI    0  2  3    C20A1 

VIN   4  0  AC  .01  ♦gain  AC  input 

VIN   4  0  PULSE (0  .01  10U  IN  IN  . 13M  .26M)  ♦slew  rate  input 
Rin   2  4  le3 
Rout  2  3  le5 

♦PLOT  ANALYSIS 

.AC  DEC  40  le3HZ  le6HZ  ♦gain  AC  analysis 

.PRINT  AC  VDB(3,0)  ♦ 

.PLOT   AC  VDB(3,0)  ♦ 

. TRAN  5U  . 3M  ♦slew  rate  analysis 

.PRINT  TRAN  V(3)  V(4)  * 

.PLOT   TRAN  V (3)   (0,-1)  V(4)   (0,.01)  * 

.END 


55 


frequency  was  determined  by  finding  the  frequency,  indexed  by  n3db  in  the  code,  whose 
data  point  came  closest  to  the  value  of  (OdB  frequency  minus  3dB).  Then  percentages  of 
gain  and  of  3dB  were  calculated  for  comparison  with  previous  results.  Since  the  op  amps 
simulated  were  different  from  those  actually  irradiated,  a  relative  or  percentage 
comparison  was  more  appropriate  and  yielded  better  results  than  a  straight  comparison. 

2.    Slew  Rate  Degradation 

As  mentioned  previously,  in  order  to  obtain  slew  rate  data,  separate  source  and 
analysis  lines  were  needed  which  were  commented  with  *  in  the  programs  presented 
below.  Again,  an  AC  analysis  was  conducted,  where  the  output  constitutes  three 
columns  of  data:  time,  Vout,  and  Vin. 

After  loading  data  using  the  same  procedure  discussed  above,  the  data  was  then 
arranged  in  a  matrix  where  rows  one  through  twenty-one  contained  single  op  amp  circuit 
data  (time,  Vout,  and  Vin)  for  each  of  seven  |3  degradation  points,  and  rows  twenty-two 
through  forty-two  contained  associated  composite  op  amp  data.  This  huge  matrix  was 
then  transferred  into  a  second  file  for  processing,  where  the  maximum  slew  rate  was 
obtained  by  taking  the  maximum  of  the  voltage  difference  divided  by  the  time  difference 
for  each  time  interval. 
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*PSPICE    CMOS    SINGLE    OPAMP    INVERTING  AMPLIFIER 
* CONFIGURATION 

. SOBCKT    CMOSOPAMP    2,1,7 
VDD    10    0    DC    5 
VSS    9    0    DC    -5 
IREF    6    9    DC    25D 

2    L=8U  W=120U 
2    L=8D   W=120U 
L=10U   W=50D 
L=10D   W=50D 
2    L=10D  W=150D 
L=10U  W=100U 
10    10    2    L=10D  W=150U 
10    10    2    L=10D  W=150D 


Ml 

4 

1 

3  10 

M2 

5 

2 

3  10 

M3 

4 

4 

9  9  1 

M4 

5 

4 

9  9  1 

M5 

3 

6 

10  10 

M6 

7 

5 

9  9  1 

M7 

7 

6 

10  10 

M8 

6 

6 

10  10 

CC 

5 

8 

10PF 

RC 

8 

7 

10K 

CL 

7 

0 

10PF 

*NO   DEGRADATION 

.MODEL    1    NMOS     (LEVEL=2    VTO=l    KP=20D    LAMBDA=0.04) 

.MODEL    2    PMOS     (LEVEL=2    VTO=-l   KP=10D    LAMBDA=0.04) 

*-.25V   DEGRADATION 

.MODEL    1    NMOS     (LEVEL=2    VTO=.  75  KP=20D   LAMBDA=0.04) 

.MODEL    2    PMOS     (LEVEL=2    VTO=-1.25  KP=10D    LAMBDA=0.04) 

*- . 5V  DEGRADATION 

.MODEL    1    NMOS     (LEVEL=2    VTO=.5  KP=20D    LAMBDA=0.04) 

.MODEL    2    PMOS     (LEVEL=2    VTO=-l . 5  KP=10U   LAMBDA=0.04) 


*-2V  DEGRADATION 

.MODEL    1    NMOS     (LEVEL=2    VTO=-2    KP=20U    LAMBDA=0.04) 
.MODEL    2    PMOS     (LEVEL=2    VTO=-3   KP=10D    LAMBDA=0.04) 
.ENDS    CMOSOPAMP 

♦INVERTING      AMPLIFIER    CIRCUIT 
XI         0    2    3         CMOSOPAMP 

VIN      4    0   AC    .01  *gain   AC    input 

VIN    4    0    PULSE (0     .01    10U    IN    IN    . 13M    .26M)     *slew   rate    input 
Rin      2    4    le3 
Rout   2    3    le5 

♦PLOT   ANALYSIS 
.AC   DEC    40    le3HZ    le6HZ  *gain   AC    analysis 

. PRINT   AC    VDB (3,0) 
.PLOT      AC    VDB (3,0) 

. TRAN    5U    . 3M  *slew   rate   analysis 

.PRINT   TRAN  V(3)    V(4) 
.PLOT      TRAN  V (3)      (0,-1)    V(4)      (0,.01) 
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*PSPICE    CMOS    COMPOSITE     (C20A1)     OPAMP    INVERTING 
♦AMPLIFIER  CONFIGURATION 

.  SDBCKT   CMOSOPAMP   2,1,7       *  refer    to   single   cxoos   op   amp  program 


.ENDS    CMOSOPAMP 

*C20A1    CIRCUIT 
.SDBCKT    C20A1    1,2,3 

XI    2    4    5    CMOSOPAMP 

X2    1    5    3    CMOSOPAMP 

Rl    4    1    IK 

R2  5  4  6. IK 
.ENDS  C20A1 


♦INVERTING   AMPLIFIER  CIRCUIT 
XI    0  2  3    C20A1 
VIN   4  0  AC  .01 

VIN  4  0  PULSE (0  .01  10U  IN  IN 
Rin   2  4  le3 
Rout  2  3  le5 


♦gain  AC  input 
13M  .26M)  *slew  rate  input 


♦PLOT  ANALYSIS 
.AC  DEC  40  le3HZ  le6HZ 
.PRINT  AC  VDB(3,0) 
.PLOT   AC  VDB(3,0) 
. TRAN  5U  . 3M 


♦gain  AC  analysis 
* 
* 
♦slew  rate  analysis 


PRINT  TRAN  V(3)  V(4) 

PLOT   TRAN  V (3)   (0,-1)  V(4)   (0,.01) 


.END 
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%MATLAB  LOADING  PROGRAM 

%THIS  PROGRAM  LOADS  DATA  OF  BIPOLAR  (BI)  JUNCTION 
% TRANSISTORS  AND  CMOS  (CM)  IN  ORDER  TO  PROCESS  GAIN  DATA 


load  b:\thesis\gain\oa  bi(:,l:2) 
load  b:\thesis\gain\oad50  bi ( : , 3 
load  b:\thesis\gain\oad90  bi ( : , 5 
load  b:\thesis\gain\oad95  bi ( : , 7 
load  b:\thesis\gain\oad96  bi ( : , 9 
load  b:\thesis\gain\oad97  bi ( : , 1 
load  b:\thesis\gain\oad975  bi ( : , 
load  b:\thesis\gain\c2oal  bi ( : , 1 
load  b:\thesis\gain\c2oald50  bi ( 
load  b:\thesis\gain\c2oald90  bi ( 
load  b:\thesis\gain\c2oald95  bi ( 
load  b:\thesis\gain\c2oald96  bi ( 
load  b:\thesis\gain\c2oald97  bi ( 
load  b:\thesis\gain\c2oald99  bi ( 

save  bi.m  bi  -ascii 


=oa; 

:4)=oad50 

:6)=oad90 

:8)=oad95 

:10)=oad96; 

l:12)=oad97; 

13:14)=oad975; 

5:16)=c2oal; 

:  ,17:18)=c2oald50 

: ,19:20)=c2oald90 

:  ,21:22)=c2oald95 

:  ,23:24)=c2oald96 

: ,25:26)=c2oald97 

:  ,27:28)=c2oald99 


load  b: 
load  b: 
load  b: 
load  b: 
load  b: 
load  b: 
load  b: 
load  b: 
load  b: 
load  b: 
load  b: 
load  b: 
load  b: 
load  b: 
load  b: 
load  b: 


\thesis\gain\cmos ;  cm(: , 1 : 2) =cmos ; 
\thesis\gain\cmosd25 ;  cm(: , 3 : 4) =cmosd25 ; 
\thesis\gain\cmosd50 ;  cm(: , 5 : 6) =cmosd50 ; 
\thesis\gain\cmosdlOO ;  cm( : , 7 : 8) =cmosdlOO ; 
\thesis\gain\cmosdl25 ;  cm(: , 9 : 10) =cmosdl25 ; 
\thesis\gain\cmosdl50 ;  cm( : ,11: 12) =cmosdl50 
\thesis\gain\cmosd200 ;  cm(: , 13 : 14) =cmosd200 
\thesis\gain\cmosd300;  cm(: , 15 : 16) =cmosd300 
\thesis\gain\cc;  cm(: , 17:18 )=cc; 
\thesis\gain\ccd25;  cm(: , 19 : 20) =ccd25 ; 
\thesis\gain\ccd50;  cm(: , 21 : 22) =ccd50 ; 
\thesis\gain\ccdlOO;  cm(: , 23 : 24) =ccdl00 
\thesis\gain\ccdl25;  cm(: ,25 : 26) =ccdl25 
\thesis\gain\ccdl50;  cm(: , 27 : 28) =ccdl50 
\thesis\gain\ccd200;  cm(: ,29 : 30) =ccd200 
\thesis\gain\ccd300;  cm(: , 31 : 32) =ccd300 


save  cm.m  cm  -ascii 
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%MATLAB  GAIN  AND  3dB  PROCESSING  PROGRAM 

%THIS  PROGRAM  FINDS  GAINS  AND  3dB  FREQUENCIES  FOR  ALL  CIRCUITS 
%WITHIN  THE  DATA  MATRIX  BI.M.  THE  CIRCUIT  WITHIN  WHICH  ALL  OP  AMPS 
%WERE  SIMULATED  IS  AN  INVERTING  AMPLIFIER  WITH  GAIN=100  USING  A 
%RESISTOR  RATIO  OF  10 OK/ IK  AND  Vin=10mV. 

%bi.m  is  a  matrix  of  n-by-m  elements  where  odd  m's  contain 
%frequency  data  for  each  circuit  and  even  m's  contain  Vout  in 
%volts. 

load  bi.m;   clear  temp; 

%DECLARE  CONSTANTS 
%Vin  in  dB  (Vin=lmV) 
Vindb=20*logl0 ( . 01) ; 

^initiate  matrix  indices  (n  x  m) 
m=l ; n=2 ; n3db=2 ; k=l ; 

%zero-pad  matrices 
gaindb=zeros (n,size (bi ,2) /2) ; 
pf3db=zeros (1 , size (bi ,2) /2) ; 
pgain=zeros ( 1 , size (bi , 2 ) /2 ) ; 
f 3db=zeros ( 1 , si  ze (bi , 2 ) /2 ) ; 
gain=zeros ( 1 , size (bi , 2 ) /2 ) ; 
gain3db=zeros ( 1 , size (bi , 2 ) /2 ) ; 
bidata=zeros (size (bi ,2) /2 , 6) ; 

%BEGIN  ITERATIONS 

%M  ITERATES  ACROSS  ALL  CIRCUITS 
for  m=l:2:27 

temp( : ,l:2)=bi ( : ,m:m+l) ; 

%CALCULATE  GAIN,  3dB  VOLTAGE  AND  3dB  FREQUENCY 
%CALCULATE  0 -HERTZ  GAIN  AS  GAIN=Vout-Vin 

gaindb ( 1 , k) =temp (1,2) -Vindb ; 
%GAIN3db  FINDS  THE  VOLTAGE  AT  WHICH  THE  3dB  FREQUENCY  OCCURS 

gain3db (k) =gaindb ( 1 , k) -3 . 0000 ; 
%N  ITERATES  DOWN  THE  DATA  WITHIN  EACH  CIRCUIT 
n3db=2 ; 

for  n=2 : (length (temp) -1) 
%CALCULATE  GAIN  FOR  EACH  FREQUENCY  INTERVAL 

gaindb (n,k) =temp (n ,2) -Vindb; 
%IF-END  LOOP  FINDS  FREQUENCY  WHICH  IS  CLOSEST 
%MATCH  TO  3dB 
if  (gain3db(k) -gaindb (n,k) )<abs (gain3db (k) -gaindb (n3db,k) ) 
n3db=n ; 
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%MATLAB  GAIN  AND  3dB  PROCESSING  PROGRAM  (contd.) 

%CALCULATE  GAIN  BANDWIDTH  PRODUCT 
% CONVERT  FROM  dB  TO  VOLTAGE 

gain(k)=10A(gaindb(l,k)/20) ; 
%RETRIEVE  3dB  FREQUENCY 

f3db(k)=temp(n3db,l)  ; 
%CALCULATE  GAIN  BANDWIDTH  PRODUCT 

gbwp(k)=f3db(k)  *gain(k)  ; 


% CALCULATE  PERCENTAGES 
%CALCULATE  PERCENTAGES  FOR  SINGLE  OAs 
if  k<8 

pf3db(k)  =  (f3db(k)  /f3db(l)  )  *100; 
pgain (k) = (gain (k) /gain (1) ) *100 ; 
%CALCULATE  PERCENTAGES  FOR  COMPOSITE  OAs 
else 

pf3db(k)  =  (f3db(k)/f3db(8)  )*100; 
pgain (k)= (gain (k) /gain (8) ) *100; 
end 

% STORE  DATA 

bidata(k,  :  )  =  [gain  (k)  pgain  (k)  gain3db  (k)  f3db  (k)pf3db  (k)  gbwp  (k)  ]  ; 

k=k+l; 

end 

save  b:bidata  bidata 
save  b : bigaindb  gaindb 
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%MATLAB  SLEWRATE  PROCESSING  PROGRAM 

%THIS  PROGRAM  GENERATES  INSTANTANEOUS  AND  MAXIMUM  SLEWRATES 

load  b:\thesis\slewrate\bis.m; 
clear  z ; 

%DECLARE  CONSTANTS 
%count 
m=l;n=2;k=0; 

%BEGIN  INTERATION 
for  m=l : 3 : 40 
z  ( :  ,l:3)=bis ( : ,m:m+2)  ; 

k=k+l ; 

for  n=2: (length (z) ) 

%DETERMINE  INSTANTANEOUS  SLEWRATES 

slewrate (k,n)=(z(n, 2) -z((n-l) ,2) ) / (z (n, 1) -z ( (n-1) ,1) ) 
%REMOVES  INFINITE  VALUES 
if  slewrate(k,n)<lelO 

posr (k ,n) =slewrate (k ,n) ; 
else  posr (k,n) =0 ; 
end 
end 

%DETERMINE  MAXIMUM  SLEWRATES 
maxsr (k) =max (abs (posr (k, : ) ) ) ;  end 
end 

save  slewrate  slewrate  -ascii 
save  maxsr  maxsr  -ascii 
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Vin.  RESULTS 

When  examining  the  results  presented  below,  it  is  important  to  keep  in  mind  that  this 
is  a  crude  simulation  in  that  it  considers  only  the  radiation  effects  on  major  transistor 
parameters;  there  are  a  multitude  of  radiation  effects  which  were  not  considered  either 
because  they  were  not  quantifiable,  or  they  were  difficult  to  analyze,  or  both. 

The  main  objective  for  both  BJT  and  CMOS  circuit  simulations  was  to  show  that 
composite  op  amps  offered  a  larger  bandwidth  than  their  single  counterparts  when 
irradiated.  For  the  BJT  circuits,  results  from  a  previous  thesis  [Ref.  1 1  ]  which  irradiated 
BJT  components  under  a  linear  accelerator  (LINAC)  offered  a  comparison.  Hence,  a 
further  objective  was  to  validate  the  previous  LINAC  results  by  comparison  with  the 
normalized  degradation  response. 

Insofar  as  CMOS  components  were  concerned,  no  previous  work  on  CMOS 
irradiation  was  available  for  comparison,  so  the  results  of  the  CMOS  simulation  are 
presented  and  offered  for  future  comparison  with  experimental  CMOS  component 
irradiation. 

A.    COMPOSITE  VERSUS  SINGLE  OP  AMPS 

1.    Gain  and  3dB  Frequency 

Figure  8.1,  which  provides  gain  degradation  curves  for  both  single  and 
composite  BJT  op  amps,  shows  the  more  ideal  response  the  composite  op  amp  offers 
over  the  single  op  amp.    In  fact,  the  composite  op  amp  bandwidth  is  larger  than  that  of 
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the  single  op  amp  by  a  factor  of  ten  for  this  application.  In  Figures  8.2  through  8.5,  the 
bandwidth  characteristics  are  broken  down  in  terms  of  3dB  frequency  and  gain 
bandwidth  product  (ft  or  GBWP).  In  these  figures,  each  op  amp  response  lies  between 
two  curves.  The  curves  are  an  outgrowth  of  (3-to-total-dose  mapping  using  Figure  4.2. 

The  CMOS  circuit  simulation  results  provided  in  Figures  8.6  through  8.10 
mirrored  those  of  the  BIT  circuit  results  discussed  above,  in  that  the  composite  op  amps 
provided  a  much  larger  bandwidth;  however,  the  bandwidth  improvement  over  the  single 
op  amps  was  slightly  smaller  at  a  factor  of  eight,  vice  the  ten  discussed  above. 

A  general  comparison  of  all  figures  depicts  the  greater  sensitivity  CMOS 
components  have  to  ionizing  radiation  compared  to  BJT  components:  the  CMOS  circuits 
withstood  a  total  dose  of  only  50k  rads(Si)  whereas  the  BJT  circuits  were  able  to 
withstand  at  least  1M  rads(Si)  and  up  to  100M  rads(Si). 

2.    Slewrates 

Figures  8.11  through  8.16  provide  slewrate  results  in  response  to  a  pulse  input 
for  all  circuits.  In  general,  both  the  BJT  and  CMOS  composite  op  amps  provided  a 
five-times-faster  slewrate  than  their  single  op  amp  counterparts,  validating  the  first 
observation  of  this  phenomenon  by  Sage,  who  conducted  the  previous  results.  This  is 
contrary  to  theory,  which  predicts  the  same  slew  rates  under  all  conditions  for  both 
composite  and  single  op  amps.  While  there  is  no  explanation  for  this  yet,  the  faster 
slewrate  has  great  implications  for  using  composite  op  amps  in  high-speed,  particularly 
digital,  applications. 
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Figure  8.1a  BJT  Single  Op  Amp  Degradation  Curves  versus  Frequency 
as  a  Function  of  Percent  Remaining  Beta 
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Figure  8.1b  BJT  Composite  Op  Amp  Degradation  Curves  versus  Frequency 
as  a  Function  of  Percent  Remaining  Beta 
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Figure  8.2  BJT  Composite  and  Single  Op  Amp  Gain  versus  Total  Dose 
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Figure  8.3  BJT  Composite  and  Single  Op  Amp  3dB  Frequency  versus  Total  Dose 
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Figure  8.4  BIT  Composite  and  Single  Op  Amp  Percent  3dB  Frequency 

versus  Total  Dose 
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Figure  8.5  BJT  Composite  and  Single  Op  Amp  Gain  Bandwidth  Product  (GBWP) 

versus  Total  Dose 
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Figure  8.6a  CMOS  Single  Op  Amp  Gain  Degradation  Curves  versus 
Frequency  as  a  Function  of  Threshold  Voltage 
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Figure  8.6b  CMOS  Composite  Op  Amp  Gain  Degradation  Curves  versus 
Frequency  as  a  Function  of  Threshold  Voltage 
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Figure  8.7  CMOS  Composite  and  Single  Op  Amp  Gain  versus  Total  Dose 
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Figure  8.8  CMOS  Composite  and  Single  Op  Amp  3dB  Frequency 

versus  Total  Dose 
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Figure  8.9  CMOS  Composite  and  Single  Op  Amp  3dB  Frequency 

versus  Total  Dose 
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Figure  8.10  CMOS  Composite  and  Single  Op  Amp  Gain  Bandwidth  Products 

(GBWP)  versus  Total  Dose 
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Figure  8.11  displays  the  maximum  slewrate  for  each  degradation  point  of  the 
BIT  circuits;  again,  the  twin  curves  are  an  outgrowth  of  P-to-total-dose  mapping  using 
Figure  4.2.  Comparing  Figures  8.12  and  8.13,  it  is  interesting  to  see  that  the  BJT  single 
op  amp  experienced  only  a  decrease  in  slewrate  for  increasing  total  dose,  whereas  the 
BJT  composite  op  amp  experienced  a  decrease  in  slewrate  accompanied  by  delays.  The 
slewrate  results  of  the  CMOS  circuits  depicted  in  Figures  8.14  through  8.16  parallel  and 
confirm  those  of  the  BJTs,  except  for  the  lower  total  dose  threshold. 

B.    COMPARISON  WITH  LINAC  RESULTS 

1.    Non  Rad  Hard  Op  Amp  Comparison 

This  simulation  was  originally  intended  only  for  comparison  with  the  non  rad 
hard  component  data  of  Figure  6.1.  The  thought  was  that,  due  to  the  simplicity  of  the 
modeled  op  amps,  the  radiation  response  of  the  basic,  inexpensive,  non  rad  hard  op  amps 
would  compare  better  than  that  of  the  complex,  costly,  rad  hard  op  amps.  Hence,  the  non 
rad  hard  data  was  primarily  the  focus  of  attention.  Referring  to  Figure  6.1,  the  relative 
gain  loss  versus  total  dose  compared  fairly  well  with  the  lower  boundary  curve  of  the 
simulation  results  in  Figure  8.2.  Both  sets  of  gains  started  decreasing  around  3M  rads(Si) 
(Figure  8.2  is  pseudo-logarithmic),  with  a  slightly  steeper  gain  degradation  rate  thereafter 
for  the  LINAC  op  amps  (Sage's  data).  The  3dB  frequency  and  slewrate  data,  however, 
did  not  compare  well.  While  the  previous  results  showed  a  marked  difference  of 
approximately  25%  in  percent  3dB  frequency  between  the  single  and  composite  op  amps, 
the  simulated  response  in  Figure  8.4  showed,  idealistically,  no  difference  in  the  two. 
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Figure  8.1 1  BJT  Composite  and  Single  Op  Amp  Slewrate  versus  Total  Dose 
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Figure  8.12  BJT  Single  Op  Amp  Slewrate  Degradation  Curves 
versus  Time  as  a  Function  of  Percent  Remaining  Beta 
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Figure  8.13  BJT  Composite  Op  Amp  Slewrate  Degradation  Curves 
versus  Time  as  a  Function  of  Percent  Remaining  Beta 
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Figure  8.14  CMOS  Composite  and  Single  Op  Amp  Slewrate  versus  Total  Dose 
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Figure  8.15  CMOS  Single  Op  Amp  Slewrate  Degradation  Curves 
versus  Time  as  a  Function  of  Percent  Remaining  Beta 
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Figure  8.16  CMOS  Composite  Op  Amp  Slewrate  Degradation  Curves 
versus  Time  as  a  Function  of  Percent  Remaining  Beta 
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Additionally,  30-60%  3dB  frequency  decrease  was  witnessed  prior  to  1M  rads(Si),  where 
the  simulated  results  revealed  no  decrease  until  approximately  3M  rads(Si). 

2.    Rad  Hard  Op  Amp  Comparison 

The  above  results  were  somewhat  depressing.  After  reinvestigating  all  models 
and  programs  to  check  for  any  error  which  would  cause  such  a  discrepancy  in  response, 
it  was  decided  that  the  reason  why  the  data  did  not  compare  very  well  was  due  to  the  fact 
that  the  modeled  op  amps  were  very  different  than  those  that  were  experimentally 
irradiated  in  the  LINAC,  and  would  be  expected  to  respond  differently.  However,  due  to 
a  sudden  and  unexpected  arabesque  analytical  leap,  it  was  determined  that  this  thinking 
was  wrong:  the  simplicity  of  the  simulation  yields  an  ideal  response  to  ionizing  radiation, 
so  the  more  ideal  the  component  which  the  simulation  is  to  be  compared  with,  or  the 
more  rad  hard  the  component,  the  closer  the  results  would  be.  Hence,  a  rad  hard  op  amp 
comparison  was  included  as  well. 

Referring  to  Figures  6.2  and  6.3  for  18M  rads(Si)  and  68M  rads(Si)  exposure, 
respectively,  the  comparison  improved  considerably  with  Figure  6.2  in  that  all  three  op 
amps  irradiated  to  18M  rads(Si)  fell  just  within  the  normalized  gain  and  3dB  frequency 
versus  total  dose  windows  in  Figures  8.2  and  8.4.  Figure  6.3  data,  where  the  op  amps 
were  irradiated  to  68M  rads(Si),  proved  to  be  an  even  better  comparison,  because  the 
data  fell  in  the  middle  of  the  above-mentioned  windows. 
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IX.  CONCLUSIONS  AND  RECOMMENDATIONS 

There  are  two  conclusions  that  may  be  gleaned  from  all  of  this.  First  and  foremost  is 
that  composite  op  amps  are  faster  and  provide  a  broader  operational  bandwidth  relative 
to  single  op  amps. 

The  second  conclusion  is  that  the  more  rad  hard  the  component  is  to  be  modeled,  the 
better  the  prediction  of  its  response  to  total  dose  radiation  will  be  as  a  result  of  this 
simulation  method.  Characteristics  which  rad  hard  devices  may  be  immune  to,  such  as 
stray  capacitances,  noise,  structure,  surface  conditions,  etc.,  which  are  difficult  to  predict, 
were  not  simulated.  However,  these  characteristics  may  drive  a  significant  portion  of  the 
non  rad  hard  component  response,  so  much  so  that  this  model  simulation  method  is  not 
very  useful  in  predicting  non  rad  hard  response  to  total  dose  radiation. 

Additionally,  it  is  hard  to  believe  that  even  rad  hard  BJT  devices  can  still  operate 
after  exposure  to  50M  rads(Si),  as  this  is  equivalent  to  50  years  in  the  densest  regions  of 
the  Van  Allen  Belts.  We  may  never  know  whether  this  is  actually  true,  since  other 
satellite  subsystems,  particularly  batteries  and  propellant,  do  not  last  that  long.  However, 
the  results  of  this  simulation  stand  behind  the  results  of  the  LINAC  irradiation. 

In  order  to  fully  understand  the  relationship  between  transistor  response  and  op  amp 
circuit  response,  a  recommendation  for  future  work  is  to  irradiate  a  chip  containing 
transistors  and  at  least  one  op  amp  made  from  those  same  transistors.  Measurements  of 
the  responses  in  these  experiments  have  to  also  be  obtained  in-situ  in  order  to  avoid 
annealing   effects.      Upon   conducting   such   an   experiment,   a   comparison   could   be 
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performed  of  the  measured  output  from  the  irradiated  op  amp  model  with  the  predicted 
response  of  the  op  amp  using  this  simulation  procedure,  but  with  new  transistor 
parameter  degradation  values  extracted  from  the  in-situ  transistor  measurements. 
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